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ABSTRACT 


The  signal-to-noise  performance  of  the  ICP  has  been 
compared  for  several  nebulization  systems  including  a 
fixed  tip  cross-flow  nebulizer,  a  concentric  nebulizer  and 
an  ultrasonic  nebulizer.  The  precision  obtainable  with 
these  nebulizers  has  been  ascertained  for  a  range  of 
integration  times  10  msec  -  30  sec.  It  was  found  that 
little  improvement  in  precision  is  obtained  with  increasing 
integration  time.  In  fact  in  the  case  of  the  ultrasonic 
nebulizer  precision  decreased  with  increasing  integration 
time. 

To  clarify  these  results  characteristic  noise  power 
spectra  were  obtained  from  the  ICP  operating  with  a  variety 
of  nebulizers.  These  noise  power  spectra  were  obtained 
by  Fourier  analysis  of  analyte  emission  signals. 

Below  5  Hz  the  noise  power  spectra  show  a  marked 
dependence  on  the  type  of  nebulizer  used.  While  a  strict 
1/f  character  was  not  observed,  the  spectra  were  clearly 
dominated  by  low  frequency  components.  Above  5  Hz  the 
noise  power  spectra  were,  broadly  speaking,  independent 
of  the  nebulizer  used.  Distinct  peaks  were  observed  in 
the  noise  power  spectra  in  the  200  Hz  to  400  Hz  region. 

The  position  and  intensity  were  dependent  on  rf  power  and 
coolant  gas  flow  rate.  In  addition,  the  intensity  of  these 
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peaks  correspond  exactly  to  similar  peaks  found  in  the 
acoustic  noise  power  spectrum  of  the  ICP .  The  cause  of 
these  peaks  is  believed  to  result  from  rotation  of  the 
plasma . 

Correlation  studies  have  been  carried  out  between 
noise  in  the  emission  signal  and  various  plasma  parameters 
such  as  incident  rf  power,  reflected  rf  power  and  nebulizer 
spray  chamber  pressure.  No  correlation  was  observed 
between  fluctuations  in  either  the  incident  or  reflected 
power  levels  and  the  emission  signal.  Correlation  was 
however,  observed  between  fluctuations  in  the  spray  chamber 
pressure  and  the  emission  signal.  These  pressure  fluc¬ 
tuations  were  caused  by  the  flow  of  liquid  droplets  in  the 
drain  hose  from  the  spray  chamber.  This  problem  was 
eliminated  by  the  use  of  a  " soak-away"  spray  chamber 
arrangement . 

The  use  of  internal  standards  for  improving  ICP 
precision  is  investigated.  Good  correlation  is  observed 
between  analyte  and  a  variety  of  internal  standard  elements. 
When  these  internal  standards  are  used, factors  of  ap¬ 
proximately  two  improvement  in  precision  are  noted. 
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CHAPTER  I 


NOISE  IN  INDUCTIVELY  COUPLED  PLASMA  EMISSION  SPECTROMETRY 

A.  Introduction  to  the  ICP 

The  inductively  coupled  plasma  (ICP)  as  a  source  for 
atomic  emission  analysis  was  first  described  by  Greenfield 
et  al.  (1)  in  1964.  A  year  later  Wendt  and  Fassel  (2) 
published  the  results  of  the  ICP  work  that  they  had  been 
concurrently  pursuing.  The  design  of  the  ICP  torch  used 
in  both  of  these  studies  was  based  on  the  earlier  design 
of  Reed  (3/4)  who  had  an  ICP  system  for  growing  refractory 
crystals . 

The  ICP  has  found  widespread  use  in  analytical 
laboratories  throughout  the  world.  The  rapidly  growing 
usage  of  the  ICP  has  spurred  the  publication  of  many 
general  review  articles,  which  may  be  consulted  for  an 
overview  of  the  technique  (5-12). 

The  factors  that  make  the  ICP  attractive  as  a  spectro¬ 
scopic  source  include  high  sensitivity,  wide  linear  dynamic 
range  and  relative  freedom  from  chemical  interferences. 
These  advantages,  coupled  with  the  ability  of  the  ICP  to 
excite  a  much  wider  variety  of  elements  than  flames,  make 
the  ICP  uniquely  suited  for  simultaneous  analysis  by 
atomic  emission  spectrometry. 
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Many  fundamental  studies  have  been  carried  out  into 
the  characteristics  of  the  ICP .  While  much  interest  has 
been  focussed  on  excitation  mechanisms  and  source  temp¬ 
erature,  comparatively  little  fundamental  work  has  been 
carried  out  on  ICP  noise.  The  purpose  of  this  work  is  to 
study  some  of  the  types,  causes  and  characteristics  of  ICP 
noise.  In  addition  a  number  of  techniques  for  improving 
ICP  precision  will  be  investigated.  However,  before 
entering  into  the  discussion  of  ICP  noise  a  general  review 
of  noise  in  atomic  emission  spectrometry  will  be  presented. 

B.  Characterization  of  Emission  Noise 

1.  Signal-to-Noise  Ratio 

Noise  may  be  defined  as  any  unwanted  fluctuation 
present  in  a  signal.  In  the  case  of  atomic  emission 
analysis  these  are  unwanted  variations  in  the  measured 
emission  intensity  of  the  source.  The  amplitude  of  the 
fluctuations  is  usually  expressed  as  either  the  r.m.s. 
noise  current  or  voltage,  or  more  often,  by  the  standard 
deviation  of  a  number  of  repeated  measurements  of  the 
signal.  If  the  number  of  measurements  is  reasonably 
large  then  the  standard  deviation  approximates  the  r.m.s. 


noise  (13)  . 


Usually,  in  atomic  spectroscopy  the  standard  deviation 
(a)  is  referenced  to  the  signal  intensity  (S)  to  give  a 
measure  of  the  relative  noise  or  signal-to-noise  ratio 
(SNR) . 


SNR  =  S/cr  (  1) 

The  other  approach  to  expressing  the  relative  noise  is  the 
percent  relative  standard  deviation  (RSD) . 

RSD  =  (O'/S)  x  100  (  2) 

Either  way,  these  two  quantities  are  useful  in  expressing 
the  precision  that  can  be  obtained  from  the  emission  source 
with  given  instrumentation  and  under  given  conditions. 

The  use  of  the  SNR  or  RSD  does  not  however  give  any 
information  as  to  the  effects  on  the  precision  if  the 
instrumentation  or  any  of  the  conditions  are  changed.  To 
characterize  completely  the  noise  in  the  system, it  is 
necessary  to  have  the  knowledge  of  its  frequency  dis¬ 
tribution.  A  representation  of  this  frequency  distribution 
is  a  plot  of  noise  power  density  as  a  function  of  frequency, 
which  is  known  as  the  noise  power  spectrum. 

2.  Noise  Power  Spectra 


As  was  previously  mentioned  the  noise  power  spectrum 


is  a  plot  of  noise  power  density  vs  frequency.  These  plots 
are  most  useful  as  they  enable  the  prediction  of  how 
instrumentation  and  measurement  parameters  will  affect 
the  precision. 

Two  methods  are  commonly  used  to  obtain  the  noise 
power  spectra  of  a  source.  The  first  makes  use  of  a 
frequency  analyzer.  In  this  approach  the  emission  signal 
obtained  by  means  of  a  suitable  detector  such  as  a  photo¬ 
multiplier  tube  is  fed  to  the  frequency  analyser. 
Essentially  this  is  a  narrow  band  tuned  amplifier  whose 
band  pass  is  swept  over  the  frequency  range  of  interest. 

The  signal  power  from  the  amplifier  is  then  plotted  as  a 
function  of  frequency  to  give  the  noise  power  spectrum. 

The  alternative  approach  is  to  use  a  Fourier  transformation 
to  convert  noise  in  the  time  domain  to  noise  power  density 
in  the  frequency  domain.  In  practice  this  is  accomplished 
by  sampling  and  recording  the  noise  (after  certain  band¬ 
width  limitations,  see  Chapter  IV)  as  a  function  of  time. 

A  discrete  Fourier  transformation  process  is  then  carried 
out  to  produce  the  noise  power  spectrum. 

With  either  of  these  techniques  the  noise  power 
spectrum  is  a  plot  of  noise  power  density,  the  units  of 
which  can  be  watts/Hz,  volts  /Hz  or  amps  /Hz.  Noise  power 
density  results,  in  the  case  of  the  frequency  analyser, 


, 


because  all  the  noise  within  the  band-pass  of  the  amplifier 
is  simultaneously  measured.  With  the  discrete  Fourier 
transformation  method  similar  considerations  must  be  made 
with  respect  to  the  resolution  of  the  transformation  (14) . 

Noise  power  spectra  will  be  extensively  utilized  in 
the  following  discussions  of  noise  types  and  in  the  main 
body  of  the  text  for  the  interpretation  of  ICP  noise. 

C.  Noise  Sources  in  Atomic  Emission  Spectrometry 

1.  Introduction 

The  various  components  that  comprise  a  typical  ICP 
atomic  emission  spectrometer  system  are  illustrated  in 
Figure  1.  While  this  configuration  is  typical  for  an  ICP 
spectrometer  most  other  emission  techniques  employ  a 
similar  arrangement  of  components.  For  this  reason  most 
of  the  following  discussion,  while  directed  principally  at 
the  ICP, will  also  be  applicable  to  other  emission  tech¬ 
niques.  All  of  the  illustrated  components  are  capable 
of  introducing  noise  into  the  system.  Because  of  this  the 
spectrometer  system  shown  in  Figure  1  will  be  briefly 
discussed  before  a  more  detailed  discussion  of  the  more 
important  noise  sources  and  types. 

In  most  light  sources  photons  are  emitted  randomly, 
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which  gives  rise  to  a  fundamental  noise  known  as  photon 
noise.  The  degree  of  photon  noise  is  related  to  the  photon 
flux  reaching  the  detector.  Because  of  this, photon  noise 
is  dependent  on  a  number  of  factors  including  the  intensity 
of  the  source,  the  collection  efficiency  of  the  imaging 
optics  and  the  throughput  of  the  monochromator.  This 
dependence  of  shot  noise  on  the  rate  of  photons  reaching 
the  detector  makes  it  more  convenient  to  include  photon 
noise  under  the  more  general  heading  of  detector  shot  noise. 

Source  flicker  noise  is  a  "non-fundamental"  noise 
present  in  the  emission  of  most  sources.  It  may  be  due  to 
a  combination  of  many  factors  including  in  the  case  of  the 
ICP,  fluctuations  in  nebuliser  efficiency,  power  output 
from  the  rf  generator,  support  gas  flow  rates  and  tur¬ 
bulences  within  the  source.  Source  flicker  is  the  dominant 
noise  source  in  the  ICP  when  working  with  concentrations 
well  above  the  detection  limits  for  most  emitting  species. 
Because  of  the  importance  of  flicker  noise  special 
attention  will  be  given  to  its  properties  and  character¬ 
istics  later  in  this  discussion. 

Systematic  noise  may  be  present  in  the  emission 
signals  of  some  spectral  sources.  These  signal  fluctua¬ 
tions  are  non-random  and  are  observed  as  a  distinct  periodic 
fluctuation  in  the  emission  intensity,  often  observed  at 


" 
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high  frequencies.  The  cause  of  these  periodic  fluctuations 
may  be  resonances  or  similar  oscillations  of  the  source 
itself. 

The  spectrometer  may  be  a  source  of  very  low  frequency 
noise  or  drifts.  Monochromators  and  polychromators  are 
subject  to  wavelength  drifts  due  to  thermal  instabilities. 
These  instabilities  may  be  minimised  by  careful  design  or 
by  thermostating  the  instrument.  They  will  not  be  discussed 
further . 

Detector  shot  noise  is  a  fundamental  noise  associated 
with  all  forms  of  spectroscopy.  Shot  noise  in  the  case  of 
photoemissive  detectors  is  a  combination  of  two  factors. 
These  are  the  random  arrival  of  photons  at  the  photo¬ 
cathode  (photon  noise)  and  the  subsequent  random  emission 
of  photoelectrons.  With  photomultipliers  additional  noise 
is  introduced  in  the  multiplication  process  of  the  dynode 
chain . 

Also  associated  with  the  detector  is  dark  current 
noise.  This  noise  is  associated  with  the  processes  that 
give  rise  to  the  dark  current  such  as  thermionic  emission. 
This  type  of  noise  often  has  many  of  the  properties  as¬ 
sociated  with  source  flicker  noise  and  is  referred  to  as 
dark  current  flicker  (.15,16). 

Another  obvious  source  of  noise  is  the  measurement 
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electronics.  However,  in  most  well  designed  spectrometer 
systems  electronics  noise  should  be  small  as  compared  to 
the  source  flicker  and  shot  noise  present.  For  completeness 
a  number  of  sources  of  electrical  noise  will  be  briefly  out¬ 
lined.  (a)  Johnson  noise  is  an  electronic  noise  associated 
with  resistors  and  is  often  encountered  with  the  large 
resistor  that  is  commonly  used  to  provide  a  voltage  output 
from  a  photomultiplier  or  phototube.  Johnson  noise  is 
caused  by  the  random  motion  of  electrons  in  resistive 
elements  resulting  from  thermal  agitation  (13) .  In  the 
case  of  photomultipliers  Johnson  noise  is  usually  in¬ 
significant  compared  to  even  pure  shot  noise  levels,  how¬ 
ever  in  the  case  of  a  phototube, due  to  the  smaller  currents 
being  measured,  Johnson  noise  can  be  very  significant. 

(b)  Other  forms  of  electrical  noise  are  associated  with 
semiconductors  and  electronic  tubes.  Due  to  the  diverseness 
and  complexity  of  the  subject  and  the  small  importance  of 
these  noises  in  the  context  of  atomic  spectroscopy  they 
will  not  be  discussed  here.  For  detailed  discussions  the 
extensive  literature  should  be  consulted  (17-21) .  (c)  The 

measurement  electronics  are  also  susceptible  to  the  pickup 
of  extraneous  electromagnetic  radiation.  These  interference 
effects  can  be  reduced,  invariably  to  within  acceptable 
levels,  by  proper  grounding  and  shielding  techniques  (13). 


Finally,  the  act  of  converting  the  signal  to  a 

numerical  value  introduces  a  noise  component  known  as 

quantising  noise.  It  has  been  shown  (22)  that  the  variance 

2 

due  to  quantising  has  a  value  of  q  /12  if  the  total 
rms  noise  present  on  the  signal  is  greater  than  q/2 ,  where 
q  is  the  quantising  interval.  If  the  rms  noise  present  is 
less  than  q/2  it  may  actually  be  desirable  to  introduce 
noise  into  the  system,  such  that  the  q/2  limit  is  just 
reached.  This  will  ensure  that  the  least  significant  bit 
of  the  analog-to-digital  converter  is  randomised.  This 
allows  signal  averaging  to  be  carried  out  which  results  in 
improvement  of  the  measurement  precision.  Digitization 
noise  and  the  technique  of  noise  addition  is  more  fully 
discussed  in  Appendix  A. 

The  more  important  noise  sources  encountered  in  atomic 
spectroscopy,  that  is  shot  noise,  source  flicker  noise  and 
systematic  source  noise,  will  now  be  discussed  in  somewhat 
more  detail. 

2 .  Shot  Noise 

The  discussion  of  shot  noise  in  this  section  will  be 
explicitly  for  photoemissive  devices  such  as  phototubes  and 
photomultiplier  tubes.  However, due  to  the  photon  noise 
component  of  shot  noise  similar  conclusions  may  be  drawn 
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for  solid  state  detectors. 

As  was  previously  pointed  out  there  are  two  cont¬ 
ributing  factors  to  shot  noise;  these  are  photon  noise  and 
photoelectron  noise.  In  the  case  of  photomultipliers  an 
additional  factor  is  also  present  caused  by  the  amplifica¬ 
tion  process  in  the  dynode  chain.  Detailed  discussions  of 
shot  noise  can  be  found  in  the  articles  by  Fried  (23) , 
Sharpe  (24)  and  in  the  general  literature  (see  for  instance 
Van  der  Ziel  (20,21)). 

In  a  phototube,  due  to  the  random  arrival  of  photons 
at  the  photocathode  and  the  subsequent  random  emission 
of  photoelectrons,  the  arrival  of  photoelectrons  at  the 
anode  is  random  and  obeys  a  Poisson  distribution  (7) . 
Therefore  the  number  of  electrons  (N)  arriving  in  time 
(t)  will  do  so  with  a  standard  deviation  equal  to  . 

This  will  result  in  an  average  signal  current  (i  )  which 

3. 

is  simply  related  to  N,  the  charge  on  the  electron  (e) 
and  the  measurement  time. 


i  =  Ne/t 
a 


3) 


The  average  root-mean-square  noise  current  (i  rms)  will 

5 

then  be  given  by: 


r  rms 
s 


4) 


From  Equations  (3)  and  (4)  the  noise  current  can  be  ex¬ 
pressed  as  a  function  of  the  signal  and  the  measurement 
period. 


i  rms 
s 


.l/2el/2/tl/2 


(  5) 


If  one  is  now  interested  not  in  the  average  noise  current 

but  in  the  standard  deviation  of  an  integrated  signal  (a  ) 

s 

the  rms  noise  current  is  integrated  over  the  entire  meas¬ 
urement  period. 


.l/2el/2tl/2 


(  6) 


The  expression  for  the  signal-to-noise  ratio  becomes: 


SNR 


il/2tl/2/el/2 


(  7) 


Equation  (7)  shows  the  important  conclusion  that  in  the 
case  of  shot  noise  it  is  possible  to  gain  a  square  root 
improvement  in  the  signal-to-noise  ratio  simply  by  extending 
the  measurement  period  (integration  time) .  Also  Equation 
(7)  reveals  a  square  root  dependence  of  the  signal-to- 
noise  ratio  on  the  signal  intensity.  Thus  it  is  possible 
to  reduce  the  effects  of  shot  noise  simply  by  increasing 
the  photon  flux  reaching  the  detector. 

The  same  conclusions  also  apply  to  the  photomultiplier 
tube  with  respect  to  the  effects  of  integration  time  and 
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photon  flux.  The  amplification  process  does  not  alter  the 
fundamental  shot  noise  component.  However,  statistical 
variations  in  the  number  of  secondary  electrons  emitted  by 
the  dynodes  for  each  impinging  primary  electron  adds  an 
additional  noise  factor.  The  overall  equation  for  the  rms 
noise  current  becomes: 


i  rms 
s 


1/2  1/2 . 1/2  ,,1/2 

ag  e  i  '  /t 


(  8) 


where  a  is  the  noise  enhancement  factor  for  the  dynode 

chain  induced  noise  and  g  is  the  gain  of  the  photomultiplier 

tube.  Typical  values  for  a  are  in  the  region  of  1.2  and 

5 

the  gain  is  in  the  order  of  10  (24) . 

The  noise  power  spectrum  of  shot  noise  is  white,  that 
is,  the  noise  power  density  P(f)  is  independent  of  fre¬ 
quency  (see  Figure  2a) .  For  a  phototube  an  estimate  of 
the  noise  power  density  P(f)  may  be  given  by: 

P(f)  =  2ei  (  9) 

ci 


For  a  photomultiplier  tube  additional  factors  for  the  gain 
and  the  dynode  noise  factor  must  be  included  (20): 


P(f)  =  2ageia 
3.  Source  Flicker  Noise 


(10) 


Source  flicker  noise  is  a  non- fundamental  noise 
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Figure  2.  Noise  power  spectra  (a)  white  noise,  (b)  1/f  no 


associated  directly  with  the  spectral  source-  It  is 
synonymously  known  as  analyte  flicker,  fluctuation  and  1/f 
noise.  The  main  characteristics  of  source  flicker  noise 
are  the  proportionality  of  the  noise  to  the  signal  and  a 
noise  power  spectrum  with  approximately  a  1/f  functionality. 

Flicker  noise  is  not  unique  to  atomic  spectroscopy. 

It  is  encountered  in  many  physical  systems  including 
electrical  noise  in  carbon  films,  semiconductors  and  zener 
diodes,  and  also  with  quartz  oscillators,  in  the  rotation 
of  the  earth  and  in  the  flow  of  highway  traffic  (17) . 

While  the  sources  of  flicker  noise  in  atomic  spectroscopy 
are  not  well  understood  the  causes  of  flicker  noise  in  many 
electronic  systems  are  well  known  and  appear  to  be 

fundamental  in  nature  and  are  associated  with  queuing 
phenomena  (17) . 

In  flame  atomic  absorption  spectrometry  Bower  and 
Ingle  (25)  have  stated  that  they  believe  source  flicker  is 
caused,  to  a  large  extent,  by  variations  in  the  rate  of 
analyte  delivery  to  the  flame  by  the  nebulizer.  As  will 
be  demonstrated  later  in  this  thesis  the  nebulizer  also 
appears  to  be  a  major  source  of  noise  in  the  ICP  system. 

Winefordner  et  al.  (26,27)  used  an  empirical  ex¬ 
pression,  that  was  previously  shown  to  be  applicable  to 
certain  electronic  tubes,  to  describe  source  flicker  noise 
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in  atomic  emission  spectrometry.  They  assumed  it  applied 
equally  well  to  evaluate  the  noise  power  density  function 
of  spectral  source  flicker  noise : 

P(f)  =  kai^/fn  (11) 

where  i^  is  the  signal,  f  is  the  frequency  and  k  is  a 
dimensionless  constant  for  a  particular  system.  This 
constant  expresses  the  ratio  of  flicker  noise  power  to 
signal  power  upon  which  the  flicker  noise  is  carried, 
a  and  n  are  also  constants  and  take  on  values  close  to  2 
and  1 . 

i)  The  Transfer  Function  and  Measurement  Effects  on 
Flicker  Noise 

It  was  possible  in  the  case  of  shot  noise  to 
calculate  on  a  statistical  basis  the  effect  of  integration 
time  on  the  signal-to-noise  ratio.  Due  to  the  non¬ 
fundamental  nature  of  source  flicker  noise  a  similar  cal¬ 
culation  is  not  possible.  However  from  a  knowledge  of  the 
flicker  noise  power  spectrum  and  the  transfer  function  of 
the  measurement  system  it  is  possible  to  predict  the  effect 
of  changing  the  integration  time. 

The  transfer  function  (H(f))  may  be  used  to  character¬ 
ize  totally  the  frequency  response  of  a  system.  It  can 


, 


be  defined  as  the  ratio  of  signal  amplitude  at  the  input 
of  a  device  (v^)  to  the  amplitude  at  the  output  (vq), 
expressed  as  a  function  of  frequency: 


H (f )  =  v  /v. 

o  i 


(12) 


in  terms  of  power  the  transfer  function  becomes: 


H  ( f ) 2  =  P  /P. 

o/  1 


(13) 


For  example  the  transfer  function  of  a  first  order  RC  low 
pass  filter  is  (13) : 


H(f)  =  1/  [  (2TrfRC)2  +  1]1/2 


(14) 


The  transfer  function  of  this  type  of  filter  is  shown  in 
Figure  3  for  a  number  of  RC  time  constants.  In  any  in¬ 
tegrated  measurement  the  transfer  function  resulting  from 
the  integration  time  itself  must  be  considered.  The 
transfer  function  associated  with  an  integration  time  (t) 
is  given  by  (13) : 


(15) 


The  effective  transfer  function  for  several  integration 
times  is  illustrated  in  Figure  4.  For  multicomponent 
systems  the  transfer  function  of  the  system  is  simply  the 
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Figure  3.  Transfer  function  of  first  order  low  pass  filter  (a) 
(b)  RC  =  0.1  s ,  (c)  RC  =  0.01  s. 
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Figure  4.  Transfer  function  of  integrator  (a) 


product  of  the  transfer  function  of  its  components: 


H  ( f )  =  H  (f )  ,  H  (f )  0  .  .  . H  ( f )  (16) 

Now  armed  with  the  knowledge  of  the  transfer 
function  of  the  measurement  system  and  the  noise  power 
spectrum  of  the  source  it  is  possible  to  calculate  the 
effect  of  integration  time  on  flicker  noise.  This  is  done 
by  combining  the  noise  power  spectrum  and  the  transfer 
function  of  the  integrator  and  integrating  over  all 
frequencies  (25) : 

i-  rms  =  [/°°P  (f)  H  (f)  2df] 1</2  (17) 

t  o 

If  one  assumes  a  noise  power  spectrum  that  is  described 
by  Winef ordner ' s  empirical  equation  (Equation  11),  setting 
a  =  2  and  n  =  1,  and  substitutes  this  expression  together 
with  the  transfer  function  of  the  integrator  into  Equation 
(17)  an  expression  for  the  average  flicker  noise  current 
(i^  rms)  can  be  derived: 

i  2  rms  =  4 (ln2)k2i  2  (18) 

r  a 

This  equation  is  more  commonly  written  in  the  form: 

if  rms  =  C  ia  (19) 

where  £  replaces  the  constant  terms  in  Equation  (18)  and 


is  known  as  the  flicker  factor.  If  the  flicker  noise  is 
now  integrated  for  a  finite  period  of  time  (t)  the  standard 
deviation  of  these  measurements  (cr^)  will  be: 

=  S  ia  t  (20) 

In  terms  of  the  signal-to-noise  ratio  the  expression 
becomes : 


SNR  =  1/e  (21) 

This  result  should  be  compared  with  the  equation  obtained 
for  shot  noise  (Equation  (7) ) .  The  important  consequences 
of  Equation  (21)  are  that  unlike  shot  noise  the  signal-to- 
noise  ratio  is  now  independent  of  the  integration  time  and 
the  photon  flux.  Thus  no  improvement  in  precision  can  be 
obtained  by  increasing  the  integration  time  or  the  signal 
intensity.  Using  this  type  of  analysis  it  is  also  possible 
to  show  that  for  a  system  in  which  the  low  frequency 
components  are  greater  than  that  predicted  by  a  1/f  re¬ 
lationship  increasing  the  integration  time  will  actually 
decrease  the  precision. 

4 .  Systematic  Source  Noise 

In  the  noise  power  spectra  of  certain  emission  sources 
distinct  peaks  are  observed.  These  peaks  are  due  to  noise 
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which  is  not  random  but  systematic  in  nature. 

Alkemade  et  al.  (28)  observed  the  presence  of  systematic 
noise  in  the  emission  of  a  premixed  H2-02~N2  flame* 

Hieftje  and  Bystroff  (29)  also  observed  similar  peaks  in 
the  noise  power  spectra  of  sheathed  air-acetylene  flames. 
Hieftje  attributed  the  origin  of  these  peaks  to  "the  swirl 
effect  of  eddy  currents  in  the  flame".  It  is  known  (30) 
that  flames  may  become  unstable  under  certain  conditions 
and  oscillate  at  high  frequencies.  Presumably,  this  could 
lead  to  similar  systematic  noise  in  flames  when  they  are 
operated  under  the  appropriate  conditions. 

Noise  of  this  type  has  been  observed  in  ICP  emission 
signals.  Walden  et  al.  (14)  observed  peaks  in  the  noise 
power  spectrum  of  the  ICP  but  was  unable  to  explain  their 
origin.  Belchamber  and  Horlick  (31)  observed  similar  peaks 
in  the  noise  power  spectrum  of  the  ICP  in  the  200  to  400  Hz 
region.  They  explained  their  origin  in  the  rotation  of  an 
assymetric  plasma  body.  The  work  leading  to  this  con¬ 
clusion  is  detailed  in  Chapter  IV. 

Because  of  the  high  frequency  of  these  systematic 
noises  it  is  possible  to  eliminate  effectively  their 
influence  by  the  use  of  extremely  short  ingetration 
times.  However,  systematic  noise  of  this  nature  may  be  a 
severe  problem  when  monitoring  transient  emission  signals  or 


* 


23 


when  using  the  plasma  as  an  atom  or  ion  cell  for  atomic 
absorption  or  fluorescence  studies.  These  noise  components 
also  pose  a  problem  when  a  Fourier  transform  spectrometer 
is  used  as  part  of  the  measurement  system  (32) . 

D.  The  Combined  Effects  of  Noise  on  the  Signal 

1.  The  Combination  of  Noises 

The  total  noise  present  at  the  output  of  the  spectro¬ 
meter  is  a  combination  of  all  noise  components  introduced 

by  the  system.  For  independent  noise  sources  the  total 
2 

variance  (a  )  is  the  sum  of  the  variances  of  the  in¬ 
dividual  noises: 


2  =  2 
°T  Q FLICK 


(22) 


+  aSHOT  +  a ELECTRON 


Care  must  be  taken  however  to  ensure  that  the  noise  sources 
are  independent,  if  they  are  not  then  an  additional  factor 
the  correlation  coefficient  (r)  must  be  introduced  (27,33): 


(23) 


The  value  of  the  correlation  coefficient  lies  in  the  in¬ 
terval  -1  to  +1  and  expresses  the  degree  of  linear  relation¬ 
ship  between  the  two  noises.  In  most  cases  the  noises 
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Equation  (22)  may  be  used  to  combine  them.  One  major  ex¬ 
ception  exists.  This  is  encountered  when  combining  flicker 
noise  in  the  emission  signal  and  flicker  noise  in  the 
background  emission.  In  this  case  it  is  reasonable  to 
expect  some  degree  of  correlation  to  exist  between  the  two 
noises.  This  calculation  may  be  necessary  under  con¬ 
ditions  of  low  signal-to-background  ratio  where  background 
fluctuations  become  the  dominant  source  of  noise  (34). 

2.  Multiplicative  and  Additive  Noise 

Noise  from  spectroscopic  emission  sources  can  either 
be  additive  or  multiplicative  in  its  effects  on  the  signal 
(29,35)  .  One  example  of  additive  noise  is  variations  in 
the  emission  intensity  of  the  background.  In  this  case  the 
noise  is  merely  superimposed  upon  the  signal.  Multiplica¬ 
tive  noise  is  a  modulation  in  the  intensity  of  the  emission 
signal.  This  modulation  may  be  brought  about  by  a  number 
of  factors  including  variations  in  the  analyte  transport 
rate  to  the  source  or  changes  in  the  excitation  conditions. 

Additive  noise  originating  in  the  background  emission 
is  obviously  more  of  a  problem  under  low  signal-to-back¬ 
ground  ratios.  Salin  and  Horlick  (34)  showed  that  in  the 
ICP  by  altering  the  rf  power  level  they  could  bring  about 
changes  in  the  signal-to-background  ratios  of  emission 
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lines  with  accompanying  changes  in  their  signal-to-noise 
ratios.  For  SrII  emission  at  a  signal-to-background  ratio 
of  22  the  signal -to -noise  ratio  was  430.  Decreasing  the 
signal-to-background  ratio  to  5  decreased  the  signal-to- 
noise  ratio  to  290  as  the  additive  background  noise  became 
more  important. 

When  operating  at  low  radiant  fluxes  the  limiting 
noise  may  be  due  to  additive  dark  current  flicker.  This 
was  shown  to  be  true  for  the  ICP  by  Greenfield  and 
Thorburn-Burns  (36) .  Under  these  conditions  the  use  of 
source  modulation  will  result  in  a  gain  in  signal-to-noise 
ratio  (37) . 

Multiplicative  noise  may  be  dealt  with  by  the  use  of 
the  internal  standard  technique.  This  technique  involves 
ratioing  the  emission  intensity  of  the  analyte  to  either 
the  emission  intensity  of  another  element  or  some  other 
parameter  whose  amplitude  is  modulated  by  similar  multi¬ 
plicative  noise.  If  the  dominant  noise  is  not  multiplica¬ 
tive  but  additive  then  it  is  obvious  that  no  gain  in 
precision  can  be  expected  by  the  use  of  an  internal 
standard,  indeed  the  precision  may  be  degraded. 

E.  Inductively  Coupled  Plasma  Noise  Studies 


While  few  fundamental  studies  have  been  carried  out 
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specifically  on  ICP  noise,  information  of  this  nature  is 
°f  critical  importance  if  the  technique  is  to  reach  its 
full  potential.  Meddings  et  al.  (38)  pointed  out  the 
great  financial  burden  often  placed  on  the  precision  of 
ICP  analysis.  While  most  studies  into  ICP  noise  have 
focussed  on  the  detection  limit  relatively  little  at¬ 
tention  has  been  paid  to  the  precision  of  measurements 
of  concentrations  considerably  in  excess  of  the  detection 
limit.  The  signal-to-background  ratio  has  been 
used  as  a  figure  of  merit  (39,40).  Salin  and  Horlick  (34) 
showed  that  this  ratio  is  only  significant  when  the 
limiting  noise  present  in  the  system  is  related  to  the 
plasma  background  emission. 

Boumans  (41,42)  stated  that  for  his  system  there 
was  a  constant  standard  deviation  of  approximately  1%  in 
his  signal  and  also  a  1%  standard  deviation  in  the  back¬ 
ground  emission  over  the  entire  spectral  range.  Later 
Boumans  (9)  concluded  that  the  increase  in  relative 
standard  deviation  with  decreasing  concentration  was  in 
keeping  with  a  system  in  which  background  carried  fluc¬ 
tuation  noise  dominates.  Greenfield  (43)  stated  emphatically 
that  noise  in  his  ICP  was  directly  proportional  to  the 
signal . 

Salin  and  Horlick  (34)  studied  the  signal-to-noise 
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ratio  performance  of  the  ICP  using  a  photodiode  array 
spectrometer.  Their  conclusions  were  that  well  above  the 
detection  limit  analyte  flicker  is  the  dominant  noise. 

As  the  concentration  is  decreased  a  gradual  transition 
takes  place  until  at  very  low  concentrations  detector 
noise  becomes  limiting.  A  similar  result  was  shown  by 
Greenfield  and  Thorburn-Burns  (36) :  under  certain  con¬ 
ditions  of  low  radiant  fluxes  the  principal  component  of 
the  background  signal  with  photomultiplier  tube  was  dark 
current . 

Ediger  et  al.  (44)  in  a  study  on  the  selection  of  ICP 
instrumental  properties  noted  that  for  arsenic,  manganese 
and  uranium  emitting  species  little  improvement  in  the 
precision  of  analysis  was  to  be  gained  by  the  use  of  long 
integration  times. 

From  these  studies  it  seems  safe  to  conclude  that  at 
concentrations  largely  in  excess  of  the  detection  limit 
the  limiting  ICP  noise  is  source  flicker.  As  the  detection 
limit  is  approached  background  flicker  noise  becomes 
dominant  and  under  some  conditions  this  background  noise 
is  not  due  to  the  plasma  but  is  detector  dark  current 
noise. 

The  causes  of  source  flicker  are  not  well  understood. 
However,  the  nebulisation  step  plays  an  important  role. 
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Greenfield  (45)  compared  ten  nebulizers  for  use  with  the 
ICP  and  found  considerable  differences  in  their  precision. 
Salin  and  Horlick  (34)  stated  that  in  their  opinion  the 
most  probable  sources  of  source  flicker  lie  in  the  nebulizer 
sample  delivery  and  gas  flow  systems  of  the  ICP.  This 
point  of  view  manifests  itself  in  the  fact  that  attempts 
to  improve  the  precision  of  the  ICP  have  largely  concen¬ 
trated  on  nebulizer  and  sample  transport  system  design. 

The  first  design  of  a  nebulizer  specifically  for  use 
with  the  ICP  was  a  variable  tip  cross-flow  nebulizer 
operating  at  low  gas  flow  rates  described  by  Knisely  et 
al.  (46)  in  1974.  Myers  et  al.  (47)  pointed  out  the  very 
considerable  effects  of  tip  alignment  with  this  type  of 
cross-flow  nebulizer,  on  the  signal-to-noise  ratio. 

Meddings  et  al.  (38,48)  concluded  that  in  terms  of  short 
term  precision  cross-flow  nebulizers  gave  the  best  results 
though  in  terms  of  long  term  stability  the  rigid  tube 
alignment  present  in  concentric  nebulizers  helped  give 
better  reproducibility.  They  then  described  a  fixed  tip 
cross-flow  nebulizer  operated  at  a  pressure  of  200  psi. 

This  nebuliser  when  used  in  conjunction  with  electronic 
mass  flow  controllers  to  control  delivery  of  argon  to  the 
ICP,  gave  a  precision  of  <0.5%  RSD  in  routine  analysis.  A 
fixed  tip  cross-flow  nebulizer  used  with  an  impact  ball 
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has  also  been  described  by  Novak  et  al.  (49)  for  use  with 
the  ICP .  The  detection  limits  obtained  with  this  nebulizer 
are  comparable  to  those  obtained  with  others.  No  data  is, 
however,  available  on  its  precision. 

Schutyser  and  Janssens  (50)  evaluated  a  number  of 
spray  chamber  designs  for  ICP  use.  These  spray  chambers 
were  modifications  to  the  original  design  of  Scott  et  al . 
(51) .  Relative  standard  deviations  using  these  spray 
chambers  with  the  same  Meinhard  concentric  nebulizer  for 
Lil  (670.8  nm)  and  Cdll  (226.5  nm)  were  in  the  range  of 
0.4%  for  the  best  case  to  2.6%  for  the  worst. 

It  is  the  intent  of  the  work  contained  in  this  thesis 
to  expand  on,  and  explore  more  fully  some  of  the  origins 
of  source  flicker  and  other  noises  associated  with  the  ICP. 
It  is  also  the  intent  to  investigate  procedures  such  as 
internal  standardisation  in  their  effectiveness  in  removing 
the  affects  of  source  flicker  noise  in  ICP  emission  analysi 
To  help  characterize  the  noise  from  the  ICP  a  study 
of  the  effect  of  integration  period  on  precision  was 
carried  out  for  a  number  of  nebulizers.  To  characterize 
the  emission  noise  more  completely  noise  power  spectra  of 
the  ICP  were  obtained,  with  a  number  of  nebulizers  and 
under  different  running  conditions. 

In  addition,  correlation  studies  have  been  carried 
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out  between  noise  in  the  emission  signal  and  fluctuations 
in  ICP  parameters  including  spray  chamber  pressure,  in¬ 
cident  power  and  reflected  power.  Correlation  studies 
have  also  been  carried  out  to  investigate  the  effectiveness 
of  a  number  of  internal  standard  elements.  Before  entering 
into  these  studies, the  ICP  system,  its  modifications  and 
the  measurement  equipment  used  to  carry  out  this  work  will 
be  described. 


CHAPTER  II 


A.  ICP  Instrumentation  and  Modifications 

1.  The  ICP  Source 

The  ICP  system  used  for  this  work  was  a  Plasma-Therm 
ICP  5000  Plasma-Therm,  Inc.,  Kresson  N.J.).  The  spec¬ 
ifications  for  this  system  are  listed  in  Table  I.  Power 
from  the  radio  frequency  generator  is  coupled  via  an 
automatic  impedance  matching  and  tuning  network  to  the  ICP 
source  (Figure  5) .  This  matching  network  allows  efficient 
power  transfer  to  the  plasma,  under  a  wide  range  of  op¬ 
erating  conditions.  A  four  turn  load  coil  surrounds  the 
top  of  the  plasma  torch  and  serves  to  couple  RF  power  to 
the  plasma  itself  (Figure  6)  . 

2.  Plasma  Torch 

Two  basic  designs  of  torch  were  used.  The  internal 
dimensions  were  close  to  those  suggested  by  Fassel  and 
Kniseley  (6).  The  difference  between  the  two  designs  were 
in  the  length  of  the  outer  quartz  tube  (Figure  7) .  The 
shorter  of  the  two  torches  was  used  for  the  major  part  of 
this  work;  the  tall  torch  being  mainly  used  for  comparison 
purposes  in  Chapter  IV. 

The  torches  consist  of  three  concentric  quartz  tubes. 
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TABLE  I 


ICP  System  Specifications 


R.F.  Generator  frequency 

27.12  MHz  (crystal 
controlled) 

Rated  output  power 

5  kW 

Envelope  ripple 

<5% 

Generator,  output  impedance 

50  ft 

Matching  network,  input  impedance 

50  ft 

Matching  network,  output  impedance 

wide  range 

Load  coil 


4  turn 
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Figure  5.  ICP  source  schematic. 
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Figure  6.  Spray  chamber  and  torch  assembly. 
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Figure  7.  (a)  Regular  torch,  (b)  Tall  torch.  All  dimen¬ 

sions  in  mm.  All  dimensions  identical  except  where  indicated. 
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The  central  tube  is  the  injector  tube  used  to  feed  the 
analyte  aerosol  stream  into  the  centre  of  the  plasma.  The 
intermediate  tube  is  the  auxiliary  gas  tube.  This  provides 
a  low  flow  rate  of  argon  which  serves  to  lift  the  plasma 
slightly  and  prevent  the  injector  tube  tip  from  melting. 

The  main  support  or  plasma  gas  is  tangentially  injected 
into  the  base  of  the  outer  tube.  This  gas  supply  serves 
not  only  to  support  the  plasma  but  is  judiciously  adjusted 
to  maintain  a  stable  plasma  discharge  and  prevent  the  outer 
quartz  tube  from  melting. 

Around  the  top  of  the  torch  and  over  the  load  coils 
a  quartz  bonnet  is  placed.  This  prevents  arcing  from  the 
plasma  body  to  the  coils  during  the  ignition  procedure. 

3.  Analyte  Delivery  System 

The  analyte  delivery  system,  essentially,  consists  of 
two  parts,  the  nebulizer  and  the  spray  chamber.  The 
nebulizer  converts  the  analyte  solution  into  a  fine  mist. 

It  is  mounted  inside  the  spray  chamber.  The  main  function 
of  the  spray  chamber  is  to  allow  via  impingement  and 
gravitational  settling  only  the  smaller  aerosol  droplets  to 
reach  the  plasma.  Typically,  less  than  10%  of  the 
nebulised  liquid  actually  reaches  the  plasma,  the  majority 
settling  out  in  the  spray  chamber.  The  spray  chamber  used 
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for  the  major  part  of  this  work  was  of  the  Scott  design 
(51)  and  is  illustrated  in  Figure  6. 

4.  Pneumatic  Nebulizers 

The  three  pneumatic  nebulizers  used  are  illustrated 
in  Figure  8.  The  glass  concentric  nebulizer  was  a  model 
TR-30-A1  supplied  by  J.E.  Meinhard  Associates,  Santa  Ana, 
California .  This  nebulizer  had  a  liquid  uptake  rate  of 
0.8  ml/min  at  a  gas  throughput  of  1  1/min,  at  a  -5  cm 
liquid  head.  The  variable  tip  cross-flow  nebulizer  was 
supplied  by  Plasma-Therm  Inc.,  Kresson,  N.J.  The  fixed 
tip  cross-flow  nebulizer  was  of  a  design  by  Blades  (52) . 
Both  designs  of  the  cross-flow  nebulizer  had  higher  fluid 
uptake  rates  than  the  concentric  nebulizer. 

All  the  pneumatic  nebulizers  were  operated,  for 
convenience,  under  conditions  of  natural  uptake. 

5.  Ultrasonic  Nebulizer 

The  ultrasonic  nebulizer,  constructed  in  this  lab¬ 
oratory,  was  of  the  batch  type  and  is  illustrated  in 
Figure  9.  It  consists  of  a  glass  tube  about  4  cm  in 
diameter  and  about  24  cm  in  length.  At  the  bottom  of  the 
tube  an  ultrasonic  transducer  is  mounted  in  a  Teflon  block. 
The  details  of  this  transducer  are  provided  in  Table  II. 
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Figure  8.  Pneumatic  nebulizers  (a)  variable  tip  cross-flow, 
(b)  fixed  tip  cross-flow,  (c)  Meinhard  concentric. 
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Figure  9.  Ultrasonic  nebulizer. 


TABLE  II 


Ultrasonic 

Transducer  Specifications 

Manufacturer 

Channel  Products  Inc. 

Model 

CPMT 

Nominal  Frequency 

1.4  MHz 

Drive  Voltage  (Max.) 

50  V  (rms) 

Drive  Current  (Max.) 

1  A  (rms) 

Dimensions 

1  1/16"  dia.  x  0.135"  thick 

Output  Mist  (Min.) 

5  cc/'min  at  2  0  W 
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About  15  ml  of  analyte  are  placed  in  the  nebulizer 
chamber  and  the  inlet  side  arm  is  closed  with  a  stopper. 
Radio  frequency  power  is  then  applied  to  the  transducer. 

A  fountain  of  analyte  rises  above  the  transducer  and  a  fine 
mist  of  analyte  droplets  is  produced.  This  mist  is  then 
swept  out  of  the  chamber  by  a  tangential  flow  of  argon 
carrier  gas.  The  carrier  gas  enters  through  a  side  arm 
about  halfway  up  the  chamber. 

The  long  path  length  from  the  transducer  to  the  exit 
of  the  chamber  and  the  baffle  also  ensure  that  to  some 
degree  the  larger  droplets  settle  out  before  the  analyte 
aerosol  leaves  the  chamber. 

The  transducer  was  powered  by  a  radio  frequency 
generator  (model  UNPS-1,  Plasma-Therm,  Inc.).  A  forward 
power  of  15  watts  was  typically  used.  The  generator  could 
be  tuned  about  its  nominal  frequency  of  1.4  MHz,  by  a  fine 
frequency  control.  This  was  adjusted  for  minimum  reflected 
power,  meaning  that  the  frequency  of  the  output  power 
matched  exactly  the  resonant  frequency  of  the  transducer. 

6.  Modifications  to  the  Gas  Supply  Controls 

Modifications  to  the  original  Plasma-Therm  gas  flow 
controls  were  made  to  enable  accurate  measurement  of  all 
the  gas  flow  rates. 


The  auxiliary  gas  supply  was  re-routed  through  the 
flowmeter  originally  provided  for  the  nebulizer  gas.  The 
nebulizer  gas  was  then  controlled  using  a  Matheson  7600 
series  flowmeter  equipped  with  a  603  tube  (0.1  to  4.8  1/min) 
The  needle  valve  on  this  flowmeter  was  positioned  at  the 
top  (exit)  of  the  rotameter.  This  arrangement  allows 
accurate  flow  rate  measurements  to  be  made  for  any  nebulizer 
back-pressure,  provided  the  pressure  of  the  gas  fed  to  the 
meter  is  kept  constant  at  its  calibration  pressure  (50 
psi) .  The  configuration  of  these  gas  flow  controls  is 
shown  in  Figure  5. 

7.  Ignition  and  Running  Conditions  of  the  ICP 

The  ICP  discharge  is  initiated  by  a  momentary  high 
voltage  spark,  generated  by  a  tesla  coil  connected  to  the 
base  of  the  torch.  The  operating  conditions  of  the  ICP 
are  listed  in  Table  III.  These  conditions  provide  a  stable 
and  reproducible  discharge,  with  ion  line  emission  in¬ 
tensity  peaking  at  about  17-18  mm  above  the  load  coil. 

B.  The  Optical  System 

The  ICP  source,  imaging  optics  and  monochromator  were 
mounted  on  an  optical  rail  system  (Figure  10) .  The  use  of 
this  rail  system  provided  precise  and  accurate  imaging  of 
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TABLE  III 

Typical  Operating  Conditions  of  the  ICP 


R.F.  Power  (Incident) 

1.5  kW 

R.F.  Power  (Reflected) 

<0.03  kW 

Plasma  gas  flow  rate 

16  1/min 

Auxiliary  gas  flow  rate 

0 . 7  1/min 

Nebulizer  gas  flow  rate 

dependent  on  nebulizer;  optimised 
for  maximum  signal  intensity 
(0.7  -  1.5  1/min) 

IMPEDANCE 
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the  plasma.  For  dual  monochromator  experiments  (Chapters 
IV  and  VII)  a  rail  bed  consisting  of  four  parallel  rails 
was  used. 

A  2:1  inverted  image  of  the  plasma  was  focussed  on  the 
entrance  slit  of  the  monochromator  by  a  10  cm  focal  length 
quartz  lens.  An  aperture  (3  mm)  was  placed  in  front  of  the 
entrance  slit,  thus  a  6  mm  slice  of  the  plasma  was  observed. 

The  monochromators  used  were  exclusively  Heath/GCA 
McPherson  EU-270  0.35  m  monochromators,  the  specifications 
of  which  are  listed  in  Table  IV.  The  detectors  were 
photomultiplier  tubes  (RCA,  type  1P28)  operated  at  600  V. 
The  high  voltage  power  supplies  for  these  photomultiplier 
tubes  were  model  244  High  Voltage  Supply,  Keithley  In¬ 
struments  Inc.  or  a  model  EU-701-30  Heath  Photomultiplier 
Module . 

Accurate  alignment  of  the  entire  optical  system  was 
achieved  by  means  of  a  helium-neon  laser  also  mounted  on 
the  optical  rail  bed. 

C.  Data  Acquisition  System 

1.  Measurement  Electronics 

The  measurement  electronics  available  are  illustrated 
in  Figure  11.  This  arrangement  of  components  was  typical 
for  most  of  the  experiments  carried  out.  Where  differences 
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TABLE  IV 

Monochromator  Specifications 


Model 

GCA  McPherson  EU  270 

Type 

Czerny-Turner 

Focal  length 

0.35  m 

Grating 

Plane  grating  1180  lines  per  mm, 
blazed  at  250  nm 

Reciprocal  linear 
dispersion 

Approximately  2.0  nm  per  mm 

Available  slit  width 

Continuously  variable  between 

5  and  2000  micrometers 

Typical  slit  width 

10  to  20  micrometers 

Available  slit  height 

12  mm 

Typical  slit  height 

3  mm  (aperture) 

3  mm  (aperture) 
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Figure  11.  Measurement  electronics. 


occur  these  will  be  pointed  out  in  the  separate  experimental 
sections . 

The  signal  from  the  photomultiplier  tube  is  first 
amplified  by  a  high  speed  current  amplifier  (model  42  7, 
Keithly  Instruments  Inc.,  Cleveland,  Ohio).  An  additional 
amplifier  (model  AM502,  Tektronix  Inc.,  Beaverton,  Oregon) 
is  used  for  two  purposes:  firstly, to  convert  the  0  to  10  V 
unipolar  output  of  the  current  to  a  ±5  V  bipolar  signal 
suitable  for  the  analog-to-digital  converter  of  the  data 
acquisition  system,  secondly,  the  amplifier  provided  ad¬ 
ditional  gain  and  band  pass  control  when  required. 

The  signal, in  addition  to  being  acquired  by  the  com¬ 
puter,  could  simultaneously  be  monitored  by  an  oscilloscope 
(model  SC502,  Tektronix  Inc.)  and  a  digital  voltmeter 
(model  DM501,  Tektronix  Inc.).  These  two  monitors  prove 
particularly  useful  when  setting  up  the  initial  plasma 
conditions  and  when  changing  the  wavelength  setting  of  the 
monochromator . 

2.  Minicomputer  System 

The  data  acquisition  system  was  based  on  a  PDP  11/10 
minicomputer  system  (Digital  Equipment  Comp.,  Maynard,  Mass.) 
This  computer  was  equipped  with  a  16k  byte  memory  core. 

This  core  size  was  sufficient  for  data  acquisition  routines 
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but  was  inadequate  for  the  larger  Fourier  transformation 
programs  used  to  produce  the  noise  power  spectra  and  cross 
correlation  functions.  These  programs  were  run  on  a  similar 
computer  system  but  having  a  larger/  32k  byte  core. 

The  data  acquisition  computer  is  equipped  with  a 
laboratory  peripheral  interface  (model  LPS  11,  Digital 
Equipment  Corp.).  This  unit  contains  a  sample  and  hold 
amplifier  and  a  12  bit  successive  approximation  analog-to- 
ditigal  converter  (ADC) .  This  ADC  system  can  be  made  to 
sequentially  read  any  of  eight  analog  input  channels.  In 
addition  this  unit  provides  the  digital  I/O  ports  that 
were  used  to  control  the  external  sample  and  hold  amp¬ 
lifiers  used  in  experiments  requiring  the  simultaneous 
monitoring  of  two  channels. 

After  signal  acquisition,  data  was  usually  stored  on 
a  1  M  byte  hard  disc  storage  system  (model  RK05,  Digital 
Equipment  Corp.).  Further  work  up  and  plotting  of  these 
data  were  carried  out  on  the  larger  computer  system. 

Plotting  of  these  data,  noise  power  spectra  and  cross 
correlation  functions  were  carried  out  on  an  incremental 
plotter  (model  130,  Zeta  Research  Inc.,  Lafayette, 

California) . 
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3. 


Dual  Channel  Data  Acquisition 


For  cross-correlation  experiments  it  is  necessary  to 
acquire  two  signals  simultaneously.  However,  the  computer 
data  acquisition  system  even  though  equipped  with  a  number 
of  analog  input  channels,  can  only  read  them  sequentially. 
This  problem  was  circumvented  by  the  use  of  two  additional 
external  sample  and  hold  amplifiers.  These  sample  and 
hold  amplifiers  (model  AD582,  Analog  Devices,  Norwood,  Mass, 
were  connected  in  parallel  (Figure  12) .  Sample  or  hold 
signals  were  provided  from  bit  0  of  the  digital  output 
port  of  the  LPS  11  unit.  This  allowed  data  to  be  acquired 
simultaneously  with  the  two  sample  and  hold  amplifiers 
and  then  to  sequentially  read  their  outputs.  The  software 
for  controlling  and  reading  the  sample  and  hold  amplifiers 
was  written  in  Macro  assembler  language  and  is  listed  in 
Appendix  C . 
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Figure  12 .  Dual  sample  and  hold  amplifiers. 


CHAPTER  III 


THE  EFFECT  OF  SAMPLE  INTEGRATION  PERIOD  ON 
MEASUREMENT  PRECISION  IN  INDUCTIVELY  COUPLED  PLASMA  EMISSION 

SPECTROMETRY  AND  RELATED  STUDIES 

A.  Introduction 

It  is  often  assumed,  almost  intuitively,  that  the 
precision  of  a  measurement  can  be  improved  by  the  increas¬ 
ing  signal  integration  period.  However,  as  was  shown  in 
Chapter  I,  the  improvement  that  can  be  achieved  is  de¬ 
pendent  on  the  exact  nature  of  the  limiting  noise  in  the 
system. 

The  signal-to-noise  ratios  will  improve  to  the  square 
root  of  the  integration  time  when  the  limiting  noise 
possesses  a  "white"  noise  power  spectrum;  such  would  be 
the  case  in  a  shot  noise  limited  system.  If,  as  in  the 
case  of  flicker  noise  domination,  the  noise  power  spectrum 
is  1/f  in  nature,  measurement  precision  cannot  be  improved 
by  increasing  the  signal  integration  period. 

The  effect  of  increasing  the  integration  period  from 
10  msec  to  30  seconds  for  ICP  emission  signals  of  several 
species  will  be  discussed.  The  signals  studied  include 
Call  (393.3  nm)  ,  Cdll  (214.4  nm)  and  Arl  (415.8  nm) 
emissions . 
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Nebulizer  type  and  design  have  proven  critical  in 
determining  the  precision  of  ICP  analysis.  The  precision 
performance  of  a  glass  concentric  (Meinhard) ,  a  variable 
tip  cross-flow  and  an  ultrasonic  nebulizer  were  evaluated 
as  a  function  of  integration  time.  To  assess  the  role  of 
the  nebulizer  in  determining  the  precision  of  ICP  emission 
signals  the  Arl  emission  signal  was  also  monitored  as  a 
function  of  integration  time.  The  emission  intensity  of 
this  line  if  the  nebulizer  gas  was  not  flowing  would  be 
independent  then  of  the  nebulization  process.  In  addition 
to  this  experiment  an  experiment  was  conducted  in  which 
mercury  vapor  was  introduced  independently  of  the  nebulizer 
into  the  carrier  gas  stream  and  the  precision  of  the 
mercury  emission  signal  measured. 

B.  Experimental 

The  Plasma-Therm  ICP  5000  system,  monochromator, 
photomultiplier  tube  and  electronic  measurement  system, 
as  described  in  Chapter  II,  was  used  for  these  experiments. 

The  amplified  signal  from  the  photomultiplier  tube 
was  digitised  and  acquired  by  the  computer  data  acquisition 
system.  The  12  bit  high  speed  ADC  (successive  approximation) 
was  clocked  at  1  kHz .  Precise  integration  periods  ranging 
from  10  msecs  to  30  seconds  were  achieved  by  summing  the 


appropriate  number  of  successive  conversions.  The 
software  for  controlling  the  integrator  was  written  in 
Macro  assembly  language  and  the  calling  program  which 
also  automatically  computed  the  signal-to-noise  ratios 
was  written  in  Fortran  IV.  These  programs  are  listed  in 
Appendix  C.  The  noise  bandwidth  of  the  system/  in 
addition  to  being  controlled  by  the  integration  time/  was 
limited  by  means  of  the  first  order  low  pass  filter  of  the 
Tektronix  amplifier;  in  all  cases  the  -3dB  point  was  set 
at  1  kHz . 

The  signal-to-noise  ratios  were  calculated  by  di¬ 
viding  the  mean  of  32  background  subtracted  signal  meas¬ 
urements  by  their  standard  deviation.  The  background 
signals  were  obtained  by  aspirating  distilled  water  and 
taking  the  mean  of  32  measurements.  The  background  was 
measured  prior  to  the  study  of  each  integration  time;  the 
intensity  was  assumed  to  remain  constant  throughout  the 
time  of  the  32  signal  intensity  measurements. 

The  results  of  these  experiments  are  presented  as 
plots  of  log  signal-to-noise  ratio  vs  log  integration 
time.  Figures  13,  14,  15,  16,  17  show  the  results  from  a 
1  ppm  solution  of  calcium  with  the  concentric  nebulizer, 
variable  tip  cross-flow  and  ultrasonic  nebulizers,  a  50  ppm 
cadmium  solution  with  the  concentric  nebulizer,  and  argon 
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emission  obtained  as  distilled  water  was  aspirated  with 
the  concentric  nebuliser. 

In  attempt  to  elucidate  the  relative  importance  of 
the  nebulisation/transport  step  in  determining  the  pre¬ 
cision  ,  an  additional  experiment  was  devised.  A  device  was 
constructed  that  allowed  mercury  vapor  to  be  efficiently 
mixed  with  the  carrier  gas  stream,  without  the  use  of  a 
nebuliser.  This  device  consisted  simply  of  a  coil  of 
glass  tubing  1  m  in  length  and  2  cm  in  diameter.  This 
was  packed  with  glass  wool  and  beads  of  metallic  mercury. 
Carrier  gas  flowed  through  the  coil  and  then  directly  via 
the  injector  tube  to  the  plasma.  It  was  necessary  to  place 
the  coil  in  a  large  beaker  of  water  for  thermal  stability. 

The  emission  intensity  of  the  Hgl  (253.7  nm)  line 
when  this  device  was  in  operation,  corresponded  closely 
to  that  obtained  when  a  100  ppm  solution  of  mercury  was 
aspirated  with  the  concentric  nebuliser.  The  precision  of 
the  results  obtained  with  mercury  vapor  introduction  device 
and  the  concentric  nebuliser  were  compared.  An  integration 
time  of  1  second  was  used.  The  signal-to-noise  ratios 
were,  again,  calculated  on  the  basis  of  32  background 
subtracted  measurements.  In  addition  each  experiment  was 
repeated,  typically  ten  times,  and  the  results  pooled  to 
provide  a  more  reliable  result.  These  results  are  re¬ 
produced  in  Table  V. 
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TABLE  V 


Comparison  of  Mercury  Vapor  Introduction  System  and  a 
Meinhard  Nebulizer  Aspirating  100  ppm  Mercury  Solution 


Mercury  Vapor  Meinhard  Nebulizer 


Signal 

0.1143 

X 

io7 

0.1107 

X 

10 

Std.  Deviation 

0.217 

X 

o 

1 — 1 

0.151 

X 

10 

Background  Signal 

0.7461 

X 

105 

0.6812 

X 

10 

Background  Std.  Deviation 

0.223 

X 

103 

0.361 

X 

10 

SNR 

543 

75.9 
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Discussion  of  Results 


The  log-log  plots  of  signal-to-noise  ratio  vs  signal 
integration  period  for  Call  (393.3  nm,  1  ppm)  is  shown  in 
Figure  13  for  the  concentric  (Meinhard)  nebuliser  and  in 
Figure  14  for  the  variable  tip  cross-flow.  For  both  the 
concentric  and  cross-flow  nebulisers  little,  if  any,  im¬ 
provement  in  signal-to-noise  ratio  is  realised  in  increasing 
the  signal  integration  period  from  10  msec  to  30  seconds 
(3.5  orders) .  This  clearly  shows  that  the  emission  signal 
from  these  two  nebulisers  possesses  a  l/fn  noise  power 
spectrum,  where  n  -  1.  Similar  results  were  obtained  for 
a  number  of  other  species,  for  instance,  the  effect  of 
integration  time  on  the  signal-to-noise  ratio  of  Cdll 
emission  (214.4  nm,  50  ppm)  is  shown  in  Figure  15.  All 
these  features  being  clear  indication  that  for  signals  well 
removed  from  the  detection  limits  the  dominant  noise  is 
source  flicker. 

In  the  case  of  the  ultrasonic  nebuliser  (Figure  16) 
the  trend  was  for  a  decrease  in  signal-to-noise  ratio  as  the 
integration  time  was  increased.  This  implies  that  a  severe 
flicker  noise  problem  exists  with  this  nebuliser  and  that 
the  low  frequency  components  of  the  noise  are  greater  than 
predicted  by  a  1/f  relationship. 


- 
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The  signal  to-noise  ratio  of  the  argon  emission  line 
(Figure  17)  was  little  affected  by  increasing  the  integra¬ 
tion  time.  The  indications  of  this  are  that  while  a  flicker 
noise  component  is  associated  with  the  nebulization/aerosol 
transportation  step  an  additional  flicker  noise  component 
is  associated  with  the  plasma  itself.  The  cause  of  this  is 
uncertain  but  may  lie  in  fluctuations  in  the  flow  rates  of 
the  support  gases. 

The  signal-to— noise  ratios  obtained  from  monitoring 
this  argon  emission  line  are  in  the  range  400  to  600, 
values  approximately  five  times  greater  than  those  obtained 
from  nebulised  species.  These  signal-to-noise  ratios 
correspond  closely  with  those  obtained  in  the  mercury  vapor 
introduction  experiment  (see  Table  V) .  These  values  of 
signal-to-noise  ratio  reveal  the  baseline  nebulizer  in¬ 
dependent  noise  of  the  ICP  system.  They  represent  the 
maximum  signal-to-noise  ratios  that  can  be  expected  from 
this  ICP  regardless  of  what  type  of  nebuliser  is  being 
used. 

D.  Conclusions 

Several  important  conclusions  can  be  drawn  from  these 
results.  Clearly,  little  is  to  be  gained  in  ICP  emission 
spectrometry  by  using  long  signal  integration  times  and  in 
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f&ct.  quite  short  integration  times  (<1  sec)  may  actually 
be  desirable.  All  these  results  strongly  imply  that  source 
flicker  noise  with  approximately  a  1/f  noise  power  spectrum 
limits  measurement  precision. 

It  is  also  apparent  that  the  majority  of  the  flicker 
noise  originates  in  the  nebuli zation/aerosol  transportation 
steps,  although  a  flicker  noise  component  is  also  associated 
with  the  plasma  itself. 

True  characterization  of  these  noise  components  re¬ 
quires  the  measurement  of  noise  power  spectra.  Such  meas¬ 
urements  have  been  carried  out  and  will  be  described  in 
Chapter  IV. 
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CHAPTER  IV 


NOISE  POWER  SPECTRA  OF  OPTICAL  AND  ACOUSTIC  EMISSION 
SIGNALS  FROM  THE  INDUCTIVELY  COUPLED  PLASMA 

A.  Introduction 

In  order  to  clarify  some  of  the  effects  noted  in 
Chapter  III, noise  power  spectra  of  the  ICP  were  measured 
under  different  operating  conditions  and  with  a  variety 
of  nebulizers. 

Noise  power  spectra  have  proven  useful  in  the  character¬ 
ization  of  noise  encountered  in  the  emission  of  a  number 
of  analytical  sources  (53-55)  and  specific  studies  towards 
both  flame  emission  (28,29,56,57)  and  atomic  absorption 
spectrometry  (25)  have  been  reported.  Noise  power  spectra 
may  allow  identification  of  noise  types,  identification  of 
noise  sources,  and,  of  course,  frequency  composition 
characterization.  Intelligent  choices  of  measurement 
conditions  such  as  signal  integration  periods  or  possible 
modulation  frequencies  may  also  be  made.  Knowledge  of 
noise  power  spectra  may  provide  sound  rationale  for  future 
system  improvement  and  the  noise  power  spectra  can  provide 
a  monitor  of  the  effectiveness  of  such  improvements. 

Finally,  noise  power  spectra  can  reveal  important  physical 
characteristics  of  a  spectral  source  that  may  otherwise 
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go  unnoticed. 

In  this  study  noise  power  spectra  of  optical  emission 
from  an  ICP  have  been  measured.  It  will  be  seen 
that  low  frequency  spectra  (0-5  Hz)  are  very  useful  for 
the  characterization  of  nebulizer  performance  and  analyte 
noise.  High  frequency  spectra  (0-500  Hz)  reveal  structure 
characteristics  of  plasma  spatial  dynamics.  Noise  power 
spectra  of  acoustic  emission  signals  and  two  channel  cor¬ 
relation  measurements  were  also  carried  out  in  order  to 
clarify  plasma  noise  characteristics.  Our  measurements 
and  observations  in  this  0-500  Hz  region  expand  on  those 
recently  reported  on  by  Walden  et  al.  (14) . 

B.  Experimental 

1.  ICP  Operating  Conditions 

The  plasma  system  was  the  Plasma-Therm  ICP  system 
detailed  in  Chapter  II.  The  operating  conditions  are 
listed  in  Table  III  unless  otherwise  stated.  All  gas  flows 
were  argon. 

2.  Data  Acquisition  and  Bandwidth  Control 

The  optical  emission  signals  after  detection  by  the 
photomultiplier  tube  were  amplified  by  means  of  a  current 
amplifier  and  an  additional  voltage  amplifier  (see  Chapter 
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II)  . 

amplification  the  bandwidth  of  the  signal  was 
limited  by  a  low  pass  filter  (model  3343 ,  Krohn-Hite 
Corporation ,  Avon,  Mass.) .  This  unit  has  two  filters  with 
8th  order  Butterworth  response.  When  these  filters  are 
cascaded  together  the  combined  effect  is  a  filter  whose 
roll-off  above  the  combined  -6  dB  point  is  -96  dB/octave. 

This  high  degree  of  bandwidth  control  is  crucial  to  ob¬ 
taining  unaliased,  and  hence,  meaningful  noise  power 
spectra . 

The  data  was  then  acquired  using  the  PDP  11/10  data 
acquisition  system.  For  the  0—500  Hz  spectra  the  sampling 
rate  was  1.4  kHz  and  the -6  dB  point  of  the  Krohn-Hite 
filter  was  set  at  700  Hz  which  is  also  the  Nyquist  fre¬ 
quency.  Without  the  use  of  this  filter,  frequencies  higher 
than  the  Nyquist  frequency  would  be  aliased  (folded) 
back  into  the  0-500  Hz  region  giving  spuriously  high  power 
information.  As  shown  in  Figure  18  the  use  of  the  filter 
reduces  aliasing  in  the  0-500  Hz  region  to  negligible 
levels.  For  the  0-5  Hz  spectra  data  were  sampled  at  a 
rate  of  14  Hz  and  the  -6  dB  point  of  the  Krohn-Hite  filter 
was  set  at  7  Hz.  The  response  curve  (transfer  function) 
of  the  filter  over  this  range  was  similar  in  shape  to  that 
for  the  0-500  Hz  region. 


Nyquist  Sampling 

Frequency  Frequency 
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3.  Computer  Calculation  of  Noise  Power  Spectra 

For  processing ,  the  data  were  transferred  to  a 

larger  PDP  11/10  computer  system  equipped  with  a  32k  byte 
core . 

The  steps  in  the  processing  of  this  data  to  obtain 
the  noise  power  spectrum  are  shown  in  Figure  19.  The  pro¬ 
grams  used  to  acquire  the  data  and  produce  the  noise  power 
spectra  are  listed  in  Appendix  C.  The  first  step  in  pro¬ 
cessing  was  d.c.  level  subtraction,  to  allow  efficient 
Fourier  transformation.  Apodization  (windowing)  of  the 
data  was  carried  out  to  prevent  the  generation  of  unwanted 
side  lobes  in  the  transformed  data.  These  side  lobes  are 
caused  by  the  convolution  of  the  rectangular  data  window 
with  the  data  itself.  Apodization  consists  of  multiplying 
the  data  with  a  mathematical  function  that  reduces  the 
effect  of  the  rectangular  data  window.  In  this  case  the 
apodization  was  accomplished  using  a  Gaussian  function 
(factor  10).  For  more  details  of  this  apodization  process 
consult  Yuen  (58).  After  apodization  the  d.c.  level  was 
again  subtracted  to  remove  any  small  level  introduced. 
Fourier  transformation  was  then  carried  out  to  convert  the 
noise  information  from  the  time  to  the  frequency  domain. 

To  accomplish  this  the  Fast  Fourier  Transform  algorithm 
(FFT)  was  used.  Finally,  signal  averaging  was  carried  out 
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Figure  19. 


Flow  chart  for  the  calculation  of  noise  power 


spectra. 


in  the  frequency  domain  to  produce  the  final  noise  power 
in  every  case  ten  individual  sets  of  transformed 
data  were  averaged  to  provide  a  more  accurate  estimate  of 
noise  power  density.  The  noise  power  spectra  were  then 
plotted.  No  attempt  was  made  to  calibrate  the  noise  power 
density  axis,  as  only  information  of  a  relative  nature  was 
required.  However,  the  units  are  proportional  to  amps2/Hz 
or  volts  /Hz . 

C.  Low  Frequency  Noise  Power  Spectra  (0-5  Hz) 

It  is  well  established  that  nebulizer  type,  adjust¬ 
ment,  and  running  conditions  can  have  a  marked  effect  on 
overall  system  performance.  The  signal  levels,  standard 
deviations,  and  signal-to-noise  ratios  for  Call  393.3 
(1  ppm)  obtained  using  four  different  nebulisers,  at  the 
time  when  the  noise  power  spectra  were  obtained,  are 
listed  in  Table  VI.  In  each  case  the  integration  time 
was  1  sec  and  the  standard  deviations  are  based  on  32 
replicate  measurements.  In  order  of  best  signal-to-noise 
ratio  performance  the  nebulizers  rank  glass  concentric 
first  followed  by  the  fixed  tip  cross-flow,  variable  tip 
cross-flow  with  the  ultrasonic  last.  However,  other  than 
this  fact,  the  signal-to-noise  ratio  provides  no  real 
guide  to  the  source  of  the  different  behaviour  or  dir- 


TABLE  VI 


Signal-to-Noise  Ratios 

Nebulizer 
Glass  Concentric 
Fixed  Tip  Cross-Flow 
Variable  Tip  Cross-Flow 
Ultrasonic 


for  Call  393.3  nm  (1  ppm) 


Signal 

a 

s/a 

14.9 

0.11 

134 

21.8 

0.25 

90 

22.9 

0.34 

72 

30.3 

1.06 

30 
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ections  for  improvement. 

Noise  power  spectra  (0—5  Hz)  for  these  same  nebulizers 
are  shown  in  Figure  20.  In  each  case  the  signal  monitored 
was  Call  393.3  (1  ppm) .  All  ordinates  have  the  same  scale. 
Interesting  and  significant  differences  between  these 
nebulizers  can  be  assessed  by  referring  to  these  noise 
power  spectra.  The  two  cross— flows  have  generally  similar 
spectra,  with  noise  power  increasing  approximately  lin¬ 
early  as  one  moves  from  5  Hz  towards  0  Hz  with  a  strong 
increase  near  0  Hz.  However,  the  fixed  tip  cross  flow 
clearly  has  generally  reduced  noise  level  at  all  frequencies 
when  compared  to  the  variable  tip  cross-flow.  The  fixed 
tip  cross-flow  was  constructed  in  this  laboratory  and 
consists  of  two  glass  capilliary  tubes  epoxied  in  a  fixed 
cross-flow  configuration  (Figure  8) .  This  provides  a 
more  stable  arrangement  than  the  adjustable  cross-flow. 

The  adjustable  cross-flow  was  the  Plasma  Therm  unit  based 
on  the  original  design  of  Kniseley  et  al.  (46)  (also 
illustrated  in  Figure  8) .  Thus  the  fixed  tip  arrangement 
for  a  cross-flow  nebulizer  seems  to  be  an  improvement  and 
several  workers  are  developing  cross-flow  nebulizers  with 
very  rigid  tip  alignment  (48,49). 

The  glass  concentric  nebulizer  (Meinhard)  has  a  very 
clean  noise  power  spectrum  with  only  a  weak  (relative  to 


Noise  Power 


Fixed  Tip  Cross-Flow 


Variable  Tip  Cross-Flqw 


Concentric 


Ultra-Sonic 


24 


a> 

5 

o 

CL 

© 

Vi 

'o 

z 


16- 


8- 


1.66  3.33 

Frequency  (Hz) 


5.00 


Figure  20.  Low  frequency  (0-5  Hz)  noise  power  spectra  of 

the  Call  393.3  nm  (1  ppm)  optical  emission  signal 
using  four  different  nebulizers. 
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the  other  nebulizers)  1/f  component  below  about  1  Hz.  The 
ultrasonic  nebulizer ,  on  the  other  hand,  has  a  complex  and 
relative  intense  noise  power  spectrum.  The  ultrasonic 
nebulizer  was  the  vertical  design  of  our  own  construction 
(Figure  9)  and  used  without  desolvation.  For  this  ultra¬ 
sonic  nebulizer  the  low  frequency  components  in  the  noise 
power  spectrum  are  intense  and  a  strong  1/f  dependence  is 
observed.  In  addition  a  peak  in  the  noise  power  spectrum 
occurs  at  about  1  Hz  and  it  is  felt  that  it  results  from 
the  rise  and  collapse  of  the  fountain  of  analyte  above  the 
transducer,  which  was  observed  to  occur  at  about  this  fre¬ 
quency.  It  is  interesting  to  note  that  the  ultrasonic 
nebulizer  actually  exhibits  lower  noise  power  levels  in 
the  2.5  to  5  Hz  region  than  either  of  the  cross  flow 
nebulisers.  If  a  unit  designed  to  control  or  eliminate 
excessive  solution  turbulence  was  employed,  its  signal-to- 
noise  ratio  performance  might  well  be  improved  over  that 
of  the  cross-flows. 

The  effect  of  analyte  concentration  on  the  low  fre¬ 
quency  noise  power  spectrum  is  illustrated  in  Figure  21. 
The  overall  intensity  of  the  noise  power  spectrum  depends 
on  analyte  concentration  (i.e., signal  level)  and  an  order 
of  magnitude  increase  in  concentration  is  accompanied  by 
approximately  a  two  order  magnitude  increase  in  noise 
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Figure_21_.  Effect  of  analyte  concentration  (Call,  393.3  nm) 
on  the  0-5  Hz  region  noise  power  spectrurr 
using  a  Meinhard  nebulizer. 


power  (i.e.  one  order  of  magnitude  in  noise  amplitude) . 
Note  also  the  rather  strong  increase  in  the  relative  im¬ 
portance  of  noise  frequencies  below  about  2  Hz  as  the 
concentration  of  analyte  increases,  all  consistant  with 
the  conclusion  that  analyte  flicker  is  the  dominant  noise. 
It  is  interesting  to  note  that  Figure  21  is  generally 
similar  to  Figure  4  of  Bower  and  Ingle  (25)  which  was 
obtained  for  atomic  absorption  measurements.  Likewise, 
their  conclusion,  based  on  the  nature  of  these  noise  power 
spectra,  that  there  is  little  merit  in  integrating  such 
signals  for  longer  than  about  10  sec  is  also  valid  for  ICP 
signals  and  in  fact,  data  presented  in  Chapter  III  il¬ 
lustrates  that  integration  periods  of  1  second  or  less 
will  normally  result  in  as  good  signal-to-noise  ratio  as 
can  be  expected. 

D.  High  Frequency  Noise  Power  Spectra  (0-500  Hz) 

For  the  most  part,  ICP  emission  signals  have  not  been 
observed  with  wide  bandwidth  measurement  systems.  With 
normal  solution  sample  introduction  using  standard  neb¬ 
ulizers  signals  are  typically  integrated  for  1  to  30  sec. 
However,  certain  ICP  measurements  require  wide  bandwidth 
measurement  systems.  In  this  laboratory  several  sample 
introduction  systems  based  on  direct  sample  insertion 
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(59) ,  laser  vaporization  (60)  and  electrothermal  vapor¬ 
ization  (61)  all  generate  transient  signal  plumes  which, 
for  basic  investigation,  must  be  time  resolved  and  hence 
observed  with  measurement  electronics  considerably  faster 
than  the  standard  integration  systems.  In  other  lab¬ 
oratories,  atomic  fluorescence  measurements  (62,63)  have 
been  carried  out  using  ICP ' s  and  the  choice  of  a  modulation 
frequency  requires  a  knowledge  of  the  ICP  noise  features 
at  higher  frequencies.  Finally,  in  this  laboratory  work 
has  been  carried  out  in  coupling  an  ICP  source  to  a  UV-VIS 
Fourier  transform  spectrometer  (64)  and  this  system  re¬ 
quires  that  the  ICP  emission  signals  be  observed  with 
measurement  system  bandwidths  ranging  from  0  to  20  kHz. 

Thus,  from  these  points  of  view  and  also  from  fundamental 
interest  a  study  was  carried  out  of  the  noise  power  spectra 
of  the  ICP  at  higher  frequencies  than  might  normally  be 
thought  of  as  important  or  of  interest. 

One  of  our  first  noise  power  spectra  in  the  region 
0-500  Hz  is  shown  in  Figure  22.  The  signal  observed  was 
Call  393.3  nm  (1  ppm).  In  addition  to  the  power  line 
features  which  are  always  present  to  varying  degrees,  a 
somewhat  broad  but  prominant  noise  feature  is  seen  at  about 
325  Hz  in  an  otherwise  relatively  flat  noise  power  spectrum. 
This  feature  was  also  observed  by  Walden  et  al.  (14)  but 
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Figure  22.  Noise  power  spectrum  (0-500  Hz)  of  the  Call 

393.3  nm  optical  emission  signal.  Meinhard 
nebulizer . 


few  details  were  presented  about  its  nature.  Preliminary 
investigation  of  the  noise  power  spectra  at  bandwidths  up 
to  2.5  kHz  (5  kHz  sampling  rate)  did  not  reveal  any  other 
interesting  features  and  efforts  were  concentrated  on 
clarification  of  the  0-500  Hz  region. 

The  effect  of  power  on  the  frequency  and  intensity  of 
this  noise  feature  is  shown  in  Figure  23.  The  frequency 
position  of  the  noise  peak  is  almost  linearly  dependent 
on  power  as  shown  in  Figure  24.  The  effect  of  coolant 
flow  rate  on  peak  frequency  is  shown  in  Figure  25. 
Increasing  the  coolant  flow  rate  from  15  to  25  L/min  at 
first  causes  the  peak  to  shift  to  higher  frequencies, 
pass  through  a  maximum  frequency  position  and  then  shift 
to  lower  frequencies  as  the  coolant  flow  rate  is  further 
increased. 

Over  the  course  of  several  weeks  of  experiments  two 
characteristic  peaks  were  often  observed  in  the  0-500  Hz 
region.  A  typical  spectrum  is  shown  in  Figure  26a.  The 
upper  peak  was  always  twice  the  frequency  of  the  lower 
peak.  Similar  noise  power  spectra  were  observed  for 
other  emitting  species  such  as  Cdll  214.6  nm  (Figure  26b) 
and  Arl  415.9  nm  (Figure  26c).  Also,  overall,  the  absolute 
and  relative  intensities  of  these  peaks  in  the  noise  power 
spectra  were  not  particularly  reproducible  from  experiment 
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Fi£ure__22.  Effect  of  plasma  power  (1  kw  -  1.9  kw)  on  the 
noise  power  spectrum  of  the  Call  393.3  nm 
optical  emission  signal.  Meinhard  nebulizer. 


Figure  24.  Plot  of  noise  feature  frequency  as  a  function 


of  power. 
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Figure  25.  Plot  of  noise  feature  frequency  as  a  function 


of  coolant  gas  flow  rate. 
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Figure  26. 


Typical  ICP  noise  power  spectra  in  the  0- 
bandwidth.  Meinhard  nebulizer. 


500  H 


to  experiment.  These  peaks  were  present  irrespective  of 
the  nebulisation  system  used.  Finally,  an  important  ob- 
servation  was  the  complete  disappearance  of  these  features 
when  a  tall  torch  (see  Figure  7)  was  used  and  emission 
observed,  at  normal  heights,  through  the  extended  tube 
(Figure  26d) . 

It  was  felt  at  this  point  that  some  aspect  of  plasma 
gas  dynamics  was  responsible  for  these  noise  power  spectrum 
features  and  characteristics.  The  most  reasonable  ex¬ 
planation  centered  on  rotation  of  a  slightly  asymmetric 
plasma  discharge  and  that  the  asymmetry  was  primarily  in¬ 
duced  by  air  entrainment  as  indicated  by  the  "tall  torch" 
experiment.  The  lower  frequency  was  that  due  to  the  fund¬ 
amental  frequency  of  plasma  rotation  (about  200  Hz  in 
Figure  26) .  Since  the  central  vertical  axis  of  the  plasma 
discharge  is  imaged  on  the  monochromator  any  asymmetry  in 
the  discharge  would  also  have  a  strong  second  harmonic 
component  as  the  discharge  rotates,  the  asymmetry  being 
viewed  twice  for  each  rotation  of  the  discharge.  This 
also  leads  one  to  conclude  that  the  absolute  and  relative 
intensities  of  these  features  in  the  noise  power  spectrum 
will  be  highly  dependent  on  the  spatial  region  of  the 
plasma  actually  observed.  Thus,  lens  focal  position, 
depth  of  discharge  viewed,  wavelength  of  observation 
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( chir omatic  aberation)  and  lateral  viewing  position  will 
affect/  to  some  extent  the  nature  of  the  noise  power 
spectrum.  This  is  shown  in  Figures  27  and  28  for  lens 
focal  position  and  lens  lateral  position.  Thus  great  care 
should  be  exercised  in  making  any  particular  conclusions 
about  the  intensities  of  these  features. 

One  other  possible  explanation  of  the  presence  of 
these  peaks  may  have  been  acoustic  resonances  in  the  torch. 
Calculations  based  on  the  well-known  organ  pipe  formula  revea 
that  the  very  lowest  frequency  that  gas  could  resonate  in 
a  torch,  of  the  dimensions  used,  would  be  about  800  Hz. 

These  calculations  agree  with  those  of  Walden  (14)  who 
also  rejected  the  possibility  that  acoustic  resonances  are 
responsible  for  the  peaks. 

E.  Acoustic  Emission  Signals 

The  plasma  discharge  is  a  good  acoustic  transducer 
and  power  modulated  plasma  discharges  have  actually  been 
used  as  loudspeakers.  Thus  it  was  felt  that  the  acoustic 
emission  signal  from  the  ICP  discharge  should  be  studied 
for  noise  features.  Using  a  simple  electret  microphone 
(Archer  270/092A)  (see  Figure  29)  the  noise  power  spectrum 
of  the  acoustic  signal  emitted  by  the  ICP  was  measured  and 
is  shown  in  Figure  30a.  Now  only  a  single  peak  is  observed 
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Figure  27 . 


spectra.  (b)  central  vertical  axis,  (a)  4  mm 
left  displacement  and  (c)  4  mm  right  displace¬ 
ment.  NO  band  head. 
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270-092A 


Figure  29. 


Microphone  circuit  used  for  acoustic  measurements. 
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in  the  noise  power  spectrum  at  a  frequency  corresponding 
to  the  lower  frequency  observed  in  the  case  of  the  optical 
emission  signals,  i.e.,at  the  frequency  of  plasma  rotation. 
This  is  further  evidence  that  the  upper  frequency  observed 
in  the  optical  emission  experiments  is  a  result  of  spatial 
observation  effects  which,  of  course,  will  not  affect  the 
acoustic  measurements.  It  is  interesting  to  note  that  the 
peak  in  the  acoustic  noise  power  spectrum  also  disappears 
when  a  tall  torch  is  utilized  (Figure  30b) .  Thus  there 
seems  to  be  a  general  acoustic  emission  signal  from  the 
ICP  at  a  frequency  corresponding  to  the  discharge  rotation 
frequency  and  caused  by  the  same  mechanism. 

F.  Simultaneous  Dual  Channel  Measurements 

The  evidence  presented  so  far  with  respect  to  the 
rotation  noise  is  still  somewhat  indirect.  An  experiment 
was  designed  that  would  allow  direct  observation  of  the 
rotation  of  the  plasma.  Two  completely  identical,  parallel 
measurement  systems  were  assembled  each  consisting  of  a 
monochromator  and  its  associated  measurement  electronics, 
see  Figure  31.  The  dual  sample  and  hold  amplifier  (Figure 
12)  was  used  to  allow  time  simultaneous  monitoring  of  the 
signals  from  both  measurement  channels  by  the  computer  data 
acquisition  system  (see  Chapter  II) .  The  optical  system 
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was  arranged  such  that  the  plasma  was  observed  simultaneously 
at  right  angles  by  the  two  monochromators  as  shown  in 
Figure  32.  High  and  low  pass  filters  were  used  to  isolate 
the  lower  frequency  of  the  two  noise  peaks  (Figure  33) 
from  the  other  frequencies  present  in  the  signal.  The 

aspirated  in  this  case  was  a  1  ppm  calcium  solution 
and  both  monochromators  monitored  the  Call  393.3  nm  emission 
signal.  Signals  from  both  channels  were  acquired  sim¬ 
ultaneously  by  the  data  acquisition  system  at  a  rate  of 
5  kHz ,  sufficient  to  oversample  the  characteristic  fre¬ 
quency  and  provide  accurate  information  about  any  relative 
phase  shifts.  The  results  shown  in  Figure  34  clearly 
indicate  that  the  signal  acquired  at  the  second  channel 
(B)  is  90°  out  of  phase  (lagging  behind)  that  acquired  by 
the  first  channel  (A) .  As  a  test,  the  signal  arising  from 
monochromator  A  was  split  at  the  output  of  the  photo¬ 
multiplier  tube  and  fed  simultaneously  through  both 
electronic  channels.  No  phase  shift  could  now  be  detected 
between  the  two  acquired  signals  (Figure  35) .  Thus  it  is 
safe  to  conclude  that  the  previously  measured  phase  shift 
is  real  and  not  an  artifact  of  the  measurement  electronics. 

If  this  measured  phase  shift  is  to  confirm  rotation  not 
only  should  the  two  signals  be  shifted  by  90°  but  the 
direction  of  the  shift  should  be  consistant  with  that  of 
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Figure  33.  Noise  peak. 
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Figure  34.  Simultaneously  digitised  noise  signal  for  Call 
393.3  nm  from  monochromator  A  (solid  line)  and 
monochromator  B  (dashed  line) ,  showing  90°  phase 
shift. 
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Figure  35.  Simultaneous  digitised  noise  signal  for  Call 
393.3  nm  from  monochromator  (A)  split  and  fed 
through  electronic  channels  (A)  and  (B) .  Show¬ 
ing  no  measurement  phase  shift. 
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the  expected  direction  of  rotation.  Signal  A  precedes 
signal  B  (Figure  34) ,  which  indicates  an  anticlockwise  dir¬ 
ection  of  rotation  that  is  consistant  with  the  direction  of 
rotation,  caused  by  the  tangential  introduction  of  plasma 
gas  (Figure  32) . 

G.  Conclusions 

The  measurement  of  noise  power  spectra  of  ICP  emission 
signals  has  allowed  the  clarification  of  several  aspects 
of  the  noise  properties  of  ICP  systems.  In  the  0  to  5  Hz 
region  noise  power  spectra  are  highly  characteristic  and 
diagnostically  useful  with  respect  to  nebulizer  performance 
and  design.  They  reveal  clearly  that  analyte  flicker  is 
the  dominant  source  of  noise. 

In  the  0-500  Hz  region  a  peak  in  the  noise  power 
spectrum  occurs  that  appears  to  result  from  rotation  of 
the  plasma  discharge  at  a  frequency  ranging  from  200  to 
400  Hz  depending  on  the  exact  experimental  conditions. 

This  feature  of  the  plasma  discharge  is  interesting  in  its 
own  right  and  is  important,  from  a  practical  point  of 
view,  when  measurements  are  carried  out  on  ICP  systems 
with  bandwidth  responses  in  this  region.  A  knowledge 
of  the  frequencies  where  these  peaks  occur  is  also  im¬ 
portant  when  source  modulation  frequencies  must  be  chosen 
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for  atomic  absorption  and  fluorescence  measurements, 
using  the  plasma  as  the  atom  or  ion  cell. 


CHAPTER  V 


CORRELATION  OF  SPRAY  CHAMBER  PRESSURE  FLUCTUATIONS 
WITH  NOISE  IN  THE  EMISSION  SIGNAL 

A.  Introduction 

In  dealing  with  source  flicker  noise  one  is  faced 
with  a  major  problem  if  it  is  necessary  to  improve  the 
precision  of  an  analysis.  Because  of  the  1/f  nature  of 
the  flicker  noise  power  spectrum,  it  is  not  possible  to 
improve  the  result  by  the  use  of  extended  integration 
periods.  Likewise,  modulation  techniques  are  of  little 
use  when  the  limiting  noise  is  source  flicker  (37)  .  Two 
alternatives  remain: 

1.  Identify  the  causes  of  source  flicker  and  either 
eliminate  or  reduce  their  effects.  This  may  be  brought 
about  by  careful  redesign  of  the  experiment  or  simply 
improvement  in  the  operation  of  one  component. 

2.  Utilize  the  internal  standard  principle.  This  requires 
the  location  of  another  parameter  whose  behaviour 
mimics  as  closely  as  possible  fluctuations  of  the 
emission  signal.  By  ratioing  the  emission  intensity 
with  this  other  parameter  the  multiplicative  effects 

of  source  flicker  noise  may  be  partially  or  totally 
removed.  Traditionally  the  other  parameter  used  for 
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internal  standardization  has  been  the  emission  in¬ 
tensity  of  another  element,  however,  there  is  no  reason 
why  any  other  physical  parameter  could  not  be  used. 

For  instance  it  is  common  practice  with  laser  vapor¬ 
isation/excitation  studies  to  ratio  the  emission 
intensity  with  the  laser  power  output. 

One  of  the  parameters  in  an  inductively  coupled 
plasma  system  that  may  be  correlatable  to  variations  in  the 
analyte  emission  intensity  is  the  pressure  in  the  nebulizer 
spray  chamber  (65) .  It  is  known  that  variations  in  the 
flow  rates  of  the  support  gases  of  the  ICP  cause  both 
intensity  and  spatial  changes  in  the  analyte  emission 
intensity  (66) .  It  is  also  to  be  expected  that  changes  in 
the  support  gas  flow  rate,  especially  the  nebulizer  gas 
will  result  in  pressure  changes  in  the  spray  chamber.  It 
seems  reasonable, bearing  these  facts  in  mind,  that  a 
correlation  may  exist  between  emission  noise  and  spray 
chamber  pressure  fluctuations. 

Preliminary  to  a  study  of  such  correlations,  a  study 
of  the  effects  of  the  three  gas  flow  rates  and  the  incident 
rf  power  on  the  spray  chamber  pressure  was  carried  out.  To  date 
few  pressure  measurements  have  been  made  associated  with  the  ICP. 
The  results  are  rather  limited  and  mainly  concern  them¬ 
selves  with  gas  flow  dynamics  in  the  plasma  itself  (67-70)  . 
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B.  Experimental  Measurement  of  Spray  Chamber  Pressure 


Pressure  measurements  of  the  interior  of  the  spray 
chamber  were  carried  out  using  an  electronic  pressure 
transducer  (model  G-605-01-2-5 ,  Gulton  Industries,  1644 
Whittier  Ave.,  Costa  Mesa,  CA  92627).  This  pressure 
transducer  was  coupled  to  the  spray  chamber  by  a  glass 
side  arm  and  a  short  length  of  Tygon  tubing  (Figure  36a) . 
Power  for  this  transducer  was  provided  by  two  6  V  batteries 
connected  in  series.  The  signal  from  the  transducer  was 
either  measured  directly  using  a  digital  voltmeter,  or 
amplified  prior  to  acquisition  by  the  minicomputer  system 
(Figure  36b) .  Prior  to  making  these  measurements  the 
pressure  transducer  was  calibrated  against  a  simple  water- 
filled  manometer. 

The  nebulizer  used  for  these  experiments  was  the 
Meinhard  concentric  glass  nebulizer.  The  Scott  type  spray 
chamber  was  used  for  the  initial  studies,  though  an  im¬ 
proved  spray  chamber  later  evolved  which  will  be  discussed 
later  in  this  chapter. 

C.  The  Effect  of  Plasma  Operating  Conditions  on  the  Spray 
Chamber  Pressure 

The  effects  of  nebulizer,  auxiliary  and  coolant  gas 
flow  rates  on  the  internal  pressure  of  the  spray  chamber 
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Figure  36. 


Pressure  transducer  for  spray  chamber  pressure 
measurement  (a)  connection  to  spray  chamber, 

(b)  electrical  connections. 


are  shown  in  Figures  37,38,39.  The  effect  of  incident 
power  is  shown  in  Figure  40.  It  should  be  noted  that  the 
pressure  axis  in  Figure  37  is  2.5  times  less  sensitive 
than  in  the  other  figures. 

When  not  a  variable  parameter,  input  power  was  kept 
constant  at  1.5  kW,  coolant  gas  flow  rate  was  16  1/min, 
auxiliary  gas  was  0.7  1/min  and  nebulizer  gas  flow  rate  was 
1.1  1/min.  All  gas  flows  were  argon.  It  should  also  be 
noted  that  the  gas  flow  controls  for  the  nebulizer  gas 
were  modified  from  the  original  Plasma-Therm  arrangement, 
to  provide  a  more  accurate  measurement  of  the  nebulizer 
gas  flow  rate  (see  Chapter  II,  in  particular  Figure  5). 

The  units  on  the  pressure  axes  of  the  plots  are  PSI  above 
atmospheric  pressure. 

As  would  be  expected,  the  most  important  parameter 
is  the  nebuliser  gas  flow  rate.  An  increase  in  the 
nebuliser  gas  flow  rate  increases  the  pressure  in  the 
spray  chamber  as  shown  in  Figure  37.  This  measurement  was 
made  with  and  without  the  plasma  lit.  The  back  pressure 
created  by  the  plasma  is  equal  to  the  difference  between 
the  two  curves.  Thus  there  is  a  back  pressure  simply  due 
to  the  passage  of  the  nebulizer  gas  through  the  injector 
tube.  An  additional  pressure  is  expected  when  the  plasma 
is  lit,  created,  perhaps  by  the  rapidly  heated  nebulizer 
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Figure  37,  Plot  of  spray  chamber  pressure  as  a  function  of 
nebulizer  gas  flow  rate. 
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Figure  38.  Plot  of  spray  chamber  pressure  as  a  function 


of  auxiliary  gas  flow  rate. 
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Figure  39.  Plot  of  spray  chamber  pressure  as  a  function 
of  plasma  gas  flow  rate. 
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Figure  40.  Plot  of  spray  chamber  pressure  as  a  function 


of  input  rf  power. 


gas  as  it  passes  through  the  injector  tube  tip  into  the 
plasma  toroid. 

Calculations  of  the  change  in  the  injector  tube 
orifice  size  as  a  function  of  temperature  (0.033%  increase 
for  a  300°C  rise)  indicate  that  this  effect  on  the  back 
pressure  is  not  important. 

In  contrast,  increased  flow  of  the  auxiliary  gas  has 
the  effect  of  decreasing  the  pressure  (Figure  38) .  This 
makes  intuitive  sense  as  an  increase  in  auxiliary  gas  flow 
rate  raises  the  plasma  body  in  the  torch  away  from  the 
injector  tube.  The  heating  of  the  injector  tube  tip  is 
thereby  lessened  as  is  the  inherent  back  pressure  of  the 
plasma. 

Coolant  (plasma)  gas  flow  rate  which  does  not  sub¬ 
stantially  alter  the  plasma  vertical  position  within  the 
torch,  has  very  little  effect  on  the  spray  chamber  pressure 
(Figure  39) . 

Finally,  increasing  the  input  power  increases  the 
spray  chamber  pressure  (Figure  40) .  This  would  be  expected 
as  the  plasma  tends  to  increase  in  length  (upwards  and 
downwards)  and  in  temperature  with  increasing  power.  This 
result  is  in  keeping  with  those  reported  by  Chase  (68) . 
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D.  Correlation  Studies 


1.  Experimental 

In  order  to  investigate  the  correlation  of  pressure 
fluctuations  in  the  spray  chamber  and  noise  in  the  emission 
signal,  it  is  necessary  to  monitor  simultaneously  the  two 
signals.  To  facilitate  this  the  signal  from  the  pressure 
transducer,  after  amplification  and  low  pass  filtering, 
was  fed  to  channel  1  of  the  dual  channel  sample  and  hold 
amplifier  (Figure  12) .  The  signal  from  the  photomultiplier 
tube,  after  amplification  and  low  pass  filtering, was  fed  to 
channel  0.  Both  signals  could  then  be  simultaneously 
acquired  by  the  sample  and  hold  amplifier  and  then  se¬ 
quentially  digitised  by  the  PDP  11/10  computer  data 
acquisition  system. 

The  bandwidth  of  both  channels  was  limited  by  the 
Krohn-Hite  low  pass  filters  (-48  dB/octave)  whose  -3  dB 
points  were  set  at  2.5  Hz.  Data  acquisition  was  at  a 
rate  of  5  Hz.  This  sampling  rate  was  appropriate  for 
observing  the  noise  characteristics  that  are  important  in 
determining  the  precision  of  analysis  with  most  low 
bandwidth  measurement  systems. 

2.  Calculation  of  the  Cross-Correlation  Function 


In  order  to  supplement  the  visual  comparison  of  the 
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two  signals  the  computer  was  used  to  calculate  their 
cross-correlation  function.  Mathematically  the  cross¬ 
correlation  function  C^CnAt)  °f  two  digitized  signals 
(a  and  b)  is  given  by  the  summation: 


C 


ab 


(nAt) 


£  a(t)  b(t  ±  nAt) 
t 


(24) 


where  At  is  the  sampling  interval  and  nAt  represents  the 
relative  displacement  (71) .  The  process  of  calculating 
this  function  directly  is  slow  in  the  computer,  due  to  the 
large  number  of  multiplications  required.  Instead  use  is 
made  of  the  relationship  (71) : 


a(t)  *  b(t)  =  Cab(T) 

I  ♦ 

Fourier  Transform  Inverse  Fourier  Transform  (25) 

I  i 

a  (f )  x  a  (t)  =  cab(f) 

It  is  quicker  to  Fourier  transform  both  digitised  signals, 
using  the  FFT  algorithm,  multiply  the  transforms  together 
and  then  carry  out  an  inverse  Fourier  transform  to  obtain 
the  cross-correlation  function.  The  scheme  used  to  do  this 
cross-correlation  was  that  described  by  Ng  and  Horlick 
(72).  A  listing  of  the  computer  programs  necessary  to 
accomplish  this  are  provided  in  reference  (73) . 
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3. 


Results  from  the  Scott  Spray  Chamber 


Figures  41  and  42  show  the  emission  signals  for  Cdll 
214.4  nm  (50  ppm)  and  Arl  415.9  nm,  fluctuations  in  the 
spray  chamber  pressure  and  the  cross-correlation  functions. 
In  this  instance  the  axis  for  pressure  fluctuations  has  not 
been  calibrated  with  respect  to  the  units  of  pressure. 

Figures  41  and  42  clearly  show  a  prominent  noise 
spike,  occurring  at  intervals  of  approximately  30  seconds 
in  the  emission  signals.  It  is  also  clear  that  a  similar, 
but  somewhat  different  in  shape,  feature  occurs  at  the  same 
time  in  the  signal  from  the  pressure  transducer.  While 
it  is  obvious  that  these  features  are  related  the  cross¬ 
correlation  function  is  of  little  extra  help  in  determining 
the  relationship. 

The  cause  of  the  noise  spikes  soon  became  apparent. 

A  long  drain  hose  (1  m)  from  the  spray  chamber  to  a  res- 
evoir  bottle  was  being  used.  Droplets  of  water,  from  the 
spray  chamber,  running  down  this  drain  hose  hit  the  sur¬ 
face  of  the  liquid  in  the  resevoir  bottle,  and  caused  the 
liquid  level  to  oscillate  slightly.  These  oscillations 
are  then  responsible  for  the  rapid  pressure  fluctuations 
which  in  turn  are  responsible  for  the  noise  spikes  in  the 
emission  signal.  By  changing  the  length  and  angle  of  the 
drain  tube  it  was  possible  to  reduce  these  effects  or  alter 
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Figure  41.  Correlation  of  Cdll  214.4  nm  (50  ppm)  emission 
with  spray  chamber  pressure  fluctuations. 
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Figure  42.  Correlation  of  Arl  415.9  nm  emission  with 


spray  chamber  pressure  fluctuations.  Scott 
spray  chamber. 
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the  repetition  rate.  It  was  not  possible,  however,  to 
eliminate  them  entirely.  A  number  of  drain  traps  were 
also  tried  but  with  somewhat  limited  success.  Finally  an 
improved  spray  chamber  design  which  eliminated  this  problem 
was  used. 

4.  "Soak  Away"  Spray  Chamber 

To  combat  the  fluctuations  in  spray  chamber  pressure 
caused  by  drain  tube  drips,  a  spray  chamber  with  an  im¬ 
proved  drain  system  was  constructed  (Figure  43) .  In  this 
design  the  Scott  type  spray  chamber  was  modified  so  that 
it  operated  in  a  vertically  upright  position.  The  con¬ 
ventional  drain  tube  was  replaced  by  a  bed  of  fine  washed 
sea  sand  supported  by  a  scintered  glass  frit.  The  layer 
of  sand  rapidly  becomes  saturated  with  analyte  providing 
a  barrier  to  the  escape  of  nebulizer  gas  but  allowing  liquid 
to  drain  rapidly  away. 

5.  Results  from  the  "Soak  Away"  Spray  Chamber 

When  the  "regular"  Scott  spray  chamber  was  replaced 
with  the  soak-away  system  the  pressure  fluctuations  dis- 
appeared  (Figures  44,45)  .  The  signal-to-noise  ratio  for  the 
measurement  of  the  Cdll  214.4  nm  (50  ppm)  signal  improved  (for 
a  1  second  integration  time)  from  73.2  to  97.9,  an  im- 
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Figure  43. 


"Soak-away"  spray  chamber. 
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Figure  44.  Correlation  of  Cdll  214.4  nm  (50  ppm)  emission 
with  spray  chamber  pressure  fluctuations. 
Soak-away  spray  chamber. 
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Figure  45.  Correlation  of  Arl  415.9  nm  emission  with  spray 


chamber  pressure  fluctuations.  Soak-away 
spray  chamber. 


provement  of  approximately  34%  over  the  original  precision. 

It  should  be  noted  now  that  no  correlation  could  be 
detected  when  this  spray  chamber  was  used  between  noise 
in  the  emission  signals  and  the  spray  chamber  pressure. 

This  was  true  for  both  the  emission  signals  of  Cdll 
(Figure  44)  and  Arl  (Figure  45) . 

E.  Conclusions  from  Spray  Chamber  Pressure  Measurement 

The  data  presented  in  Figures  37-40  indicate  that 
three  main  plasma  parameters  nebulizer  gas  flow  rate, 
auxiliary  gas  flow  rate  and  power,  significantly  effect 
the  pressure  in  the  spray  chamber-  The  other  parameter 
studied,  coolant  gas  flow  rate,  had  a  much  lesser  effect. 

The  measurement  and  correlation  of  spray  chamber 
pressure  and  emission  intensity  fluctuations  have  been  use¬ 
ful  in  pinpointing  a  problem  associated  with  the  drain  hose 
arrangement.  This  led  in  turn  to  the  development  of  an 
improved  spray  chamber  design  that  eliminated  the  need 
for  a  drain  hose. 

When  the  effect  of  droplets  running  down  the  drain 
hose  was  removed  no  measurable  correlation  could  then  be 
detected  between  spray  chamber  pressure  fluctuations  and 
noise  in  the  emission  signals. 
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CHAPTER  VI 


CORRELATION  OF  NOISE  IN  THE  INCIDENT  AND  REFLECTED  POWER 
CIRCUITS  WITH  NOISE  IN  THE  EMISSION  SIGNAL 

A.  Introduction 

A  parameter  of  critical  importance  in  influencing  the 
precision  of  ICP  analysis  is  the  stability  of  the  rf  power 
supply.  Kniseley  et  al.  (74)  reported  severe  modulation  of 
emission  line  intensities  caused  by  the  use  of  an  in¬ 
adequately  filtered  power  supply.  Ohls  (75)  suggested 
that  the  feed  back  of  Arl  emission  intensity  may  be  used 
to  stabilise  the  rf  power  output  of  the  generator. 

In  addition  to  variations  in  the  incident  power,  small 
changes  in  reflected  power  levels  can  often  be  observed. 

One  cause  of  these  small  changes  appears  to  be  variations 
in  the  rate  of  introduction  of  analyte  into  the  plasma. 

It  was  decided  to  monitor  fluctuations  in  both  the 
incident  and  reflected  power  levels  and  to  correlate,  if 
possible, these  fluctuations  with  those  in  the  emission 
intensity . 

B.  Experimental  Monitoring  of  Incident  and  Reflected 
Power  Levels 

The  incident  and  reflected  power  levels  were  monitored 
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by  the  use  of  the  circuit  illustrated  in  Figure  46.  The 
inputs  of  the  incident  or  reflected  power  meters  on  the 
front  panel  of  the  matching  and  tuning  network  were  dis¬ 
connected.  These  meters  are  current  meters  directly  con¬ 
nected  to  the  power  coupling  system  and  provide  a  current 
path  to  ground.  Any  attempt  to  block  this  path  to  ground 
adversely  effects  the  impedance  matching  characteristics 
of  the  system.  Therefore  an  operational  amplifier  (model 
25,  Princeton  Applied  Research  Corp.,  Princeton,  N.J.)  was 
used  in  the  low  impedance  current  amplifier  mode  to  convert 
this  current  to  a  voltage  signal.  This  voltage  produced 
was  proportional  to  the  square  root  of  the  power.  An 
analog  multiplier  (model  230,  Princeton  Applied  Research 
Corp.)  was  then  used  to  square  this  signal. 

Signals  from  the  analog  multiplier  were  then  low  pass 
filtered  and  acquired  by  the  minicomputer  data  acquisition 
system. 

Cross-correlation  studies  (as  described  in  Chapter  V) 
were  carried  out  between  fluctuations  in  both  the  incident 
and  reflected  power  circuits  and  in  the  emission  signal. 

In  addition  the  noise  power  spectra  of  both  the  incident 
and  reflected  power  signals  were  obtained. 
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Figure  46.  Circuit  for  monitoring  incident  and  reflected 


power  fluctuations. 
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C.  Discussion  of  Results 

The  noise  power  spectra  of  noise  in  the  incident  and 
reflected  power  circuits  are  shown  in  Figures  47  and  48. 
Simultaneously  acquired  signals  from  Call  (393.3  nm,  1  ppm) 
and  Arl  (415.9  nm)  emission  and  the  incident  and  reflected 
power  circuits  together  with  their  corresponding  correlation 
functions  are  shown  in  Figures  49-52.  The  presence  of  60 
Hz  envelope  ripple  on  the  27  MHz  output  of  the  generator 
is  clearly  shown  in  Figure  47.  Also  several  of  the  higher 
harmonics  of  60  Hz,  which  are  occasionally  visible  in 
emission  noise  power  spectra,  are  present.  These  features 
are  also  observed  in  the  reflected  power  noise  power 
spectrum.  It  seems  reasonable  to  explain  the  presence  of 
the  60  Hz  peak  present  in  the  emission  noise  power  spectrum 
as  resulting  from  the  60  Hz  modulation  of  the  plasma  by 
the  RF  generator.  It  is  also  worthwhile  to  note  that  none 
of  the  other  distinct  features  noted  in  the  emission  noise 
power  spectrum  are  seen  in  either  the  incident  or  reflected 
power  noise  spectra. 

The  cross-correlation  studies  show  that  no  correlation 
(in  the  0  -  2.5  Hz  range)  exists  between  emission  noise 
and  incident  power  fluctuations  for  either  Call  (Figure  49) 
or  Arl  (Figure  50) .  Correlation  is  also  absent  between 
emission  intensity  of  Call  (Figure  51) ,  Arl  (Figure  52) 
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Figure  49. 


Correlation  of  Call  emission  signal  with 
incident  power. 
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Correlation  of  Arl  emission  signal  with 
incident  power. 
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Figure  51. 


Correlation  of  Call  emission  signal  with 
reflected  power. 
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Figure  52. 


Correlation  of  Arl  emission  signal  with 
reflected  power. 
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and  reflected  power  fluctuations. 

From  this  study  it  can  be  clearly  seen  that  stability 
of  this  RF  generator  would  not  be  a  problem  with  most  low 
bandwidth  measurement  systems.  With  high  bandwidth 
measurement  systems  (as  discussed  in  Chapter  III)  the 
presence  of  60  Hz  and  its  harmonics  may  cause  a  decrease 
in  precision. 

Due  to  the  lack  of  any  detectable  (low  bandwidth) 
correlation  between  either  the  incident  or  reflected 
power  fluctuations  and  emission  intensity  noise  no  use 
could  be  made  of  either  of  these  parameters  as  an  internal 
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CHAPTER  VII 


THE  USE  OF  INTERNAL  STANDARDS  WITH  THE  ICP 
A.  The  Internal  Standard  Principle 

The  internal  standard  principle  was  introduced  in 
1925  by  Gerlach  (76) .  The  technique  involves  ratioing 
the  emission  intensity  of  the  analyte  to  that  of  another 
element  which  is  present  in  a  fixed  or  known  concentration. 
This  procedure  of  using  the  ratio  instead  of  the  actual 
emission  intensity  to  calculate  the  concentration  has  a 
number  of  potential  advantages.  Fluctuations  of  the 
excitation  conditions  within  the  source,  ideally,  should 
be  compensated  for  by  proportional  fluctuations  both  in 
the  analyte  and  internal  standard  emission  intensity. 

The  internal  standard  also  compensates  for  variations  in 
the  quantity  of  sample  being  presented  for  analysis  or  for 
variations  in  the  rate  of  sample  introduction  to  the  source. 

In  order  to  ensure  a  similar  behaviour  between  the 
analyte  and  the  internal  standard  a  number  of  rules  have 
been  built  up  over  the  years  for  the  choice  of  a  good 
internal  standard  element.  Ahrens  and  Taylor  (77)  sum¬ 
marise  these  rules,  which  for  convenience  are  reproduced 
here. 

1.  If  the  internal  standard  is  to  be  added  to  the  sample. 
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2. 


its  original  concentration  should  be  negligibly  low. 

If  added  to  the  sample,  the  internal  standard  should 
be  in  a  high  state  of  purity  with  respect  to  the 
elements  being  determined. 

3.  The  rates  at  which  internal  standard  and  analysis 
element  volatilize  should  be  very  similar. 

4 .  The  atomic  weights  of  the  internal  standard  and  the 
analysis  element  should  not  differ  greatly  if  either 
is  a  light  element. 

5.  The  ionization  energies  of  the  internal  standard  and 
analysis  elements  should  be  roughly  the  same. 

6.  Both  the  internal  standard  and  analysis  lines  should 
have  the  same  excitation  energies. 

7.  Both  lines  should  be  free  of  self  absorption. 

8.  When  photographic  recording  methods  are  used  both 
lines  should  be  approximately  the  same  wavelength. 

9.  The  internal  standard  line  should  if  possible  behave 
in  such  a  way  that  the  intensity  ratio  of  the  analysis 
pair  is  comparitively  insensitive  to  changes  of  matrix 
and  composition. 

Barnett,  Fassel  and  Kniseley  (78)  discussed  quantitatively 

the  importance  of  matching  some  of  these  factors.  In 

addition  they  also  include  the  partition  function,  which 

is  temperature  dependent  and  should  be  considered  when 
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trying  to  match  the  internal  standard  to  the  analyte. 

These  rules,  however,  have  been  built  up  around  the 
classical  techniques  of  emission  analysis  such  as  the  dc 
arc  or  the  ac  spark.  When  the  ICP  is  used  for  excitation 
it  may  be  possible  to  relax  some  of  these  criteria.  Self 
absorption  is  rarely  a  problem  with  the  ICP  and  photo¬ 
graphic  recording  is  now  seldom  used  for  quantitative 
purposes.  Also,  because  it  is  apparent  that  the  major 
source  of  ICP  flicker  noise  is  the  nebulization  step  it 
may  be  possible  to  relax  the  rules  for  matching  ionization 
energies,  excitation  energies,  partition  functions  and 
atomic  weights. 

Internal  Standards  have  met  with  some  success  when 
used  with  the  ICP.  Benton-Jones  (79)  used  molybdenum  as 
an  internal  standard  in  the  analysis  of  plant  tissue  ash. 

He  stated  that  the  choice  of  molybdenum  was  more  for 
convenience  than  by  design.  Watters  and  Norris  (80) 
observed  that  the  emission  intensity  variation  of  the  ICP 
of  1  or  2%  relative  standard  deviation  could  be  reduced 
to  0.5%,  by  concurrent  ratioing  of  the  intensities  to  the 
proper  internal  standard  element.  Salin  and  Horlick  (34) 
investigated  the  use  of  SrII  and  Arl  emission  as  internal 
standards  for  calcium.  They  found  that  the  signal-to-noise 
ratio  improved  by  a  factor  of  approximately  two  when 
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strontium  was  used  but  decreased  by  a  factor  of  two  when 
the  argon  line  was  used.  Meddings  et  al.  (38)  used 
scandium  (1  ppm)  as  an  internal  standard  for  the  analysis 
of  twenty  five  elements.  They  reported  an  improvement  in 
the  relative  standard  deviation  of  their  analysis  from 
about  1.1%  to  0.5%  when  a  variable  tip  cross-flow  nebulizer 
was  used.  However,  internal  standardization  brought  about 
only  an  improvement  from  0.3%  RSD  to  0.2%  RSD  when  their 
fixed  tip  cross-flow  nebulizer  was  used.  Uchida  et  al. 

(81)  used  yttrium  as  an  internal  standard  with  a  micro¬ 
sampling  technique  to  compensate  for  changes  in  the  amount 
of  sample  fed  to  the  plasma. 

B.  Some  Statistics  Associated  with  the  Use  of  an  Internal 
Standard 

When  an  internal  standard  is  used,  the  calculation  of 
concentration  (c)  is  based  on  the  ratio  of  the  emission 
intensity  of  the  analyte  (x)  to  that  of  the  internal 
standard  (y)  and  a  constant  of  proportionality  (k)  (78) . 

c  =  kx/y  (26) 

The  noise  in  the  ratio  of  the  emission  signals  (a^,  )  is 

related  to  the  noise  in  the  analyte  and  internal  standard 

signals  (a  )  and  (a  ) ,  the  average  intensity  of  the  signals 
x  y 
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and  the  correlation  coefficient  between  the  signals 
(r)  (33). 
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The  correlation  coefficient  expresses  a  measure  of  the 
linear  relationship  between  fluctuations  in  the  analyte 
and  internal  standard  signals.  The  value  of  the  correlation 
coefficient  lies  in  the  interval  -1  to  +1.  A  value  of  +1 
for  the  correlation  coefficient  would  indicate  a  perfect 
match  of  noise  in  the  internal  standard  signal  with  noise 
in  the  analyte  signal;  under  these  conditions  the  right 
hand  side  of  Equation  (27)  reduces  to  zero,  which  indicates 
that  the  internal  standard  has  successfully  eliminated 
the  noise. 

The  correlation  coefficient  may  be  calculated  from 
the  covariance  of  the  two  signals  (v  )  and  their  two 
standard  deviations  (33) . 
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The  covariance  is  calculated: 
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By  making  the  assumption  a  /x  -  a  /y  (not  too  un- 

x  y 

reasonable  in  the  case  of  the  ICP)  Equation  (27)  then 
becomes : 
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If  an  internal  standard  brings  about  an  improvement  in  the 
precision  of  an  analysis  then: 
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If  this  criterion  is  substituted  into  Equation  (31;: 
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then 


r  >  0.5  (34) 

This  shows  that  the  correlation  coefficient  between  the 
analyte  and  internal  standard  signals  must  be  greater  than 
0.5  before  the  internal  standard  brings  about  any  im¬ 
provement  in  precision.  If  the  correlation  coefficient  is 


less  than  0.5  then  the  use  of  the  internal  standard  will 
result  in  a  decrease  in  the  precision. 

Equation  (27)  is  also  useful  in  predicting  the  effect 
of  shot  noise  on  the  precision  when  an  internal  standard 
is  being  used.  At  best,  the  internal  standard  can  only 
be  expected  to  remove  the  effects  of  flicker  noise.  There 
may  be  a  high  degree  of  correlation  between  the  flicker 
noise  present  in  the  analyte  and  internal  standard  channel 
but  other  noises  such  as  shot  noise  are  of  necessity, 
independent  and  therefore  uncorreiated .  Thus,  in  con¬ 
ditions  where  flicker  noise  is  dominant  the  signals  may  be 
well  correlated  (r  +  1)  and  the  internal  standard  will  be 
very  effective  in  improving  the  precision.  However,  as 
other  noises  (such  as  shot  noise)  become  dominant  the 
correlation  coefficient  decreases  in  value,  until  it  is 
less  than  0.5  when  the  use  of  the  internal  standard  will 
degrade  the  precision. 

The  other  case  where  the  internal  standard  will  be  of 
little  use  is  under  conditions  of  low  signal-to-background 
ratio  (discussed  in  Chapter  I) .  Under  these  conditions 
fluctuations  in  the  background  emission  intensity  will 
have  an  additive  effect  on  either  the  analyte  or  internal 
standard  signals.  These  additive  fluctuations  will  result 
in  an  apparent  decrease  in  the  correlation  coefficient 
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between  the  two  signals.  For  this  reason  it  is  important 
when  using  internal  standards  to  ensure  high  signal-to- 
background  ratios  for  both  the  analyte  and  internal  standard 
emission  lines. 

C.  Experimental 

To  carry  out  the  investigation  into  the  use  of  in¬ 
ternal  standards  the  dual  monochromator  arrangement  (Chapter 
IV)  was  utilised.  In  this  system  two  identical  mono- 
chromators  were  mounted  at  90°  to  one  another,  to  view 
independently  the  same  region  of  the  plasma  (Figure  32) . 

This  arrangement  allowed  complete  flexibility  in  the  choice 
of  the  analyte  and  internal  standard  lines  that  were  to  be 
observed.  Data  from  both  monochromators  were  then  acquired 
simultaneously  using  the  dual  channel  measurement  system 
illustrated  in  Figure  31. 

The  data  acquisition  rate  was  5  Hz  with  the  -3  dB 
point  of  the  Krohn-Hite  low  pass  filters  (-48  dB/octave 
^"oil- off)  set  at  2.5  Hz.  Normally  1024  data  points  were 
collected  from  which  the  cross-correlation  function  and 
the  correlation  coefficient  were  calculated.  In  addition 
signal-to-noise  ratios  were  calculated  for  both  the 
analyte  emission  signals  and  the  ratio  of  the  analyte 
emission  intensity  to  that  of  the  internal  standard. 
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The  emission  lines  picked  for  this  study  were  selected 
from  the  tables  of  both  Boumans  and  Bosveld  (40),  and 
Winge,  Peterson  and  Fassel  (82) .  These  lines  were  sel¬ 
ected  on  the  basis  of  their  intensity  and  that  no  spectral 
overlap  problems  were  likely.  The  lines  used  and  the 
concentrations  of  the  solutions  used  are  listed  in  Table 
VII . 

D.  Results  From  the  Use  of  an  Internal  Standard 

The  signals  acquired  simultaneously  from  a  number  of 
analyte  and  internal  standard  elements  and  the  correspond¬ 
ing  cross-correlation  functions  are  shown  in  Figures  53-60. 
The  correlation  coefficient,  the  signal-to-noise  ratio  of 
the  analyte  emission  ratioed  to  that  of  the  internal  standard 
and  the  improvement  in  the  signal-to-noise  ratio  brought 
about  by  the  use  of  the  internal  standard  are  listed  in 
Table  VIII. 

When  the  same  emission  line  is  used  for  both  the 
analyte  and  the  internal  standard  the  precision  is  seen  to 
improve.  The  correlation  coefficients  for  Call/Call 
(Figure  53),  Cal/Cal  (Figure  58)  and  Mnll/Mnll  (not  plotted) 
were  0.9241,  0.8351  and  0.8082.  These  were  a  little 
disappointing  as  a  much  higher  degree  of  correlation  was 
expected,  considering  the  same  emission  line  was  being 
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TABLE  VII 


Spectral  Lines  and  Concentrations  Used  for 
Internal  Standard  Study 


nm 

ppm 

Call 

393.36 

1 

Cal 

422.67 

100 

Mnll 

257.61 

100 

Arl 

415.85 

- 

SrII 

407.77 

2 

Ini 

451.13 

110 

Sell 

361.38 

20 

Amplitude  x  108  A  Emission  Intensity  A  Emission  Intensity 
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Figure  53.  Correlation  of  Call  emission  signals  from  dual 


monochromator  system. 
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Figure  54.  Correlation  of  Call  and  Arl  emission  signals. 
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Figure  55.  Correlation  of  Call  and  SrII  emission  signals. 
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Figure  56.  Correlation  of  Call  and  Sell  emission  signals. 
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Figure  57.  Correlation  of  Call  and  Ini  emission  signals. 
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Figure  58.  Correlation  of  Cal  emission  signals  from  dual 


monochromator  systems. 
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Figure  59.  Correlation  of  Cal  and  Arl  emission  signals. 
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Figure  60.  Correlation  of  Cal  and  Ini  emission  signals. 
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monitored  simultaneously  by  the  two  monochromators.  The 
concentrations  used  were  high  enough  that  additive  back¬ 
ground  noise  should  not  be  a  problem  nor  should  shot  noise. 
The  only  remaining  factor  that  is  presently  perceived  as 
having  an  influence  on  these  results  may  be  the  90° 
viewing  of  the  plasma.  This,  however,  may  not  be  the  full 
story.  Salin  and  Horlick  (34)  reported  similar  results 
from  the  study  of  two  Call  lines  (393.3  and  396.8  nm) . 

In  the  ICP  these  lines  behave  almost  identically.  When 
the  signal-to-noise  ratio  obtained  from  one  of  these 
lines  was  compared  to  that  obtained  from  the  ratioed 
intensities,  only  a  factor  of  approximately  two  improvement 
was  noted  for  the  ratioed  intensities.  These  results  were 
obtained  using  a  photodiode  array  spectrometer,  which 
necessitates  viewing  the  plasma  from  the  same  angle  for 
both  the  analyte  and  the  internal  standard.  Therefore,  it 
seems  likely  that  some  other  factor  may  have  a  part  to 
play  in  limiting  the  ultimate  performance  of  the  internal 
standard . 

Other  elements  when  used  as  internal  standards  gen¬ 
erally  provided  modest  gains  in  the  signal-to-noise  ratio 
performance.  Scandium  and  indium,  two  elements  widely 
used  as  internal  standards  were  compared.  Strontium  was 
also  investigated  as  an  internal  standard  for  calcium 
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because  Call  and  SrII  have  remarkably  similar  spatial 
emission  characteristics  (66) . 

The  cross-correlation  functions  reveal  reasonable 
correlation  of  Call  emission  signals  with  those  of  SrII, 
Sell  and  Ini  (Figures  55,  56  and  57) .  Correlation  co¬ 
efficients  between  these  signals  had  values  of  0.9055, 
0.9165  and  0.6790.  This  led  to  modest  improvements  in  the 
signal-to-noise  ratio  by  a  factor  of  1.93  and  3.01  when 
SrII  and  Sell  were  used  as  internal  standards.  The  rel¬ 
atively  poorer  correlation  between  Call  emission  and  Ini 
emission  led  to  only  an  improvement  of  factor  of  1.03 
(or  essentially  no  change  in  the  signal  to  noise  ratio) 
when  it  was  used  as  an  internal  standard.  The  poor  per¬ 
formance  of  Ini  as  an  internal  standard  for  Call  may  be 
due  to  the  quite  different  behaviour  of  neutral  and  ionic 
species  in  the  plasma  (66) .  However,  a  better  correlation 
was  observed  between  Cal  emission  noise  and  Sell  noise 
than  between  Cal  and  Ini  (Figure  60) .  On  a  basis  of  the 
difference  in  behaviour  between  neutral  and  ionic  species 
the  reverse  affect  would  be  predicted.  Scandium  also 
proved  to  be  a  better  internal  standard  for  Mnll  emission 
than  indium.  The  correlation  coefficients  of  the  Mnll 
signal  with  Sell  and  Ini  signals  were  0.8818  and  0.5676. 
This  resulted  in  a  slight  improvement  in  the  signal-to- 
noise  ratio  by  a  factor  of  1.17  when  scandium  was  used  as 
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the  internal  standard  but  when  indium  was  used  the  signal- 
to-noise  ratio  was  degraded  by  about  twenty  percent. 

The  emission  intensity  of  Arl  was  also  monitored  in 
an  effort  to  assess  its  possible  use  as  an  internal  stand¬ 
ard.  Also,  there  has  been  some  interest  in  monitoring 
argon  emission  intensity  and  using  feed  back  from  this 
signal  to  control  rf  power  input  to  the  plasma  (75) . 

The  simultaneously  acquired  Call  and  Arl  signals,  and  the 
Cal  and  Arl  signals  are  shown  in  Figures  54  and  59.  The 
sharp  negative  peaks  present  in  both  cross-correlation 
functions  indicate  a  marked  degree  of  inverse  correlation 
between  the  analyte  and  the  argon  signals.  Negative  cor¬ 
relation  coefficients  were  also  calculated  (see  Table  VIII)  . 
When  argon  emission  intensity  was  ratioed  with  that  of  the 
analyte,  as  would  be  expected  from  the  correlation  co¬ 
efficient,  a  decrease  in  precision  was  noted  for  Call, 

Cal  and  Mnll .  These  results  are  in  keeping  with  those  of 
Salin  and  Horlick  (34)  who  also  observed  a  degree  of  neg¬ 
ative  correlation  between  Call  and  Arl  emission  noise. 

This  negative  correlation  may  result  from  variations  in 
the  nebuliser  flow  rate.  As  the  rate  increases  more  aerosol 
is  transported  to  the  plasma  so  that  the  analyte  emission 
intensity  increases.  The  other  effect  of  the  in¬ 
creased  aerosol  gas  flow  rate  is  to  blow  a  larger  "hole" 
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through  the  centre  of  the  plasma  which  results  in  a  de¬ 
crease  in  the  argon  emission  intensity.  Thus  it  is  pos¬ 
sible  to  visualise  a  likely  cause  of  the  inverse  correlation 
of  noise  in  the  analyte  and  argon  emission  signals. 

E.  Conclusions 

From  the  results  presented  it  is  possible  to  eliminate 
some  of  the  effects  of  source  flicker  noise  by  the  use 
of  a  suitable  internal  standard.  Generally,  the  gain 
in  precision  are  modest  (a  factor  of  two  being  typical) . 

While  it  seems  advisable  to  choose  for  ionic  emission 
lines  an  ion  line  as  an  internal  standard,  other  factors 
may  have  to  be  considered.  The  results  presented  here 
are  limited  to  only  a  few  elements  and  emission  lines; 
internal  standardisation  does  seem  an  area  suitable  for 
more  extensive  study. 
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CHAPTER  VIII 


SUMMARY 

From  the  experimental  results  obtained,  it  is  quite 
clear  that  the  dominant  emission  noise  of  the  ICP  under 
normal  operating  conditions  is  source  flicker.  Experiment¬ 
ally  a  major  source  of  this  flicker  noise  was  shown  to  be 
the  nebulization  step. 

The  effect  of  integration  time  on  the  precision  of 
ICP  analysis  was  investigated.  No  major  change  in  pre¬ 
cision  was  noted  for  an  increase  in  integration  time  from 
0.01  to  30  seconds.  This  was  true  for  a  number  of  neb¬ 
ulization  systems  and  for  an  ultrasonic  nebulizer  the 
precision  actually  worsened  with  increasing  integration 
time.  These  results  are  in  keeping  with  the  approximately 
1/f  nature  of  the  noise  power  spectrum  associated  with 
source  flicker. 

Noise  power  spectra  were  obtained  from  the  ICP 
operating  with  a  variety  of  nebulizers.  While  approximately 
a  1/f  noise  power  spectrum  was  observed  for  all  the  neb¬ 
ulizers  the  exact  nature  of  the  noise  power  spectrum  in 
the  0-5  Hz  region  was  dependent  on  the  nebulizer  used.  In 
the  0-500  Hz  noise  power  spectra  sharp  peaks  are  present 
at  60  Hz,  caused  by  ripple  in  the  rf  generator  output. 
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In  addition  rather  broad  but  well  defined  peaks  are 
located  in  the  200-400  Hz  range.  The  exact  frequency  that 
these  peaks  occur  at  is  dependent  on  the  rf  power  and 
plasma  gas  flow  rate.  The  origin  of  these  peaks  appears 
to  be  in  the  rotation  of  the  plasma.  This  type  of  noise 
while  not  being  a  problem  with  low  bandwidth  integrating 
type  measurement  systems  may  be  a  problem  with  wider 
bandwidth  systems.  In  techniques  requiring  source  modu¬ 
lation  such  as  atomic  fluorescence  and  absorption,  care 
should  be  taken  to  avoid  choosing  a  modulation  frequency 
close  to  these  peaks.  This  high  frequency  noise  is  also  a 
problem  when  dealing  with  high  speed  transient  signals 
and  when  Fourier  transform  spectrometers  are  used. 

Correlation  studies  between  emission  signals  and 
incident  and  reflected  power  levels  indicate  that  the  rf 
generator  (under  normal  conditions)  is  not  a  major  source 
of  ICP  noise,  apart  from  introducing  some  60  Hz  ripple. 
Reflected  power  level  fluctuations  do  not  bear  a  meas¬ 
urable  correlation  to  noise  in  the  emission  signal. 

Monitoring  the  spray  chamber  pressure  revealed  sig¬ 
nificant  pressure  fluctuations  caused  by  improper  design 
of  the  spray  chamber  drain  arrangement.  These  pressure 
fluctuations  correlated  extremely  well  with  similar 
fluctuations  in  the  emission  signals.  A  "soak-away" 


drain  arrangement  is  described  which  alleviates  the  spray 
chamber  drain  problem  and  increases  the  precision  of 
analytical  results. 

The  use  of  internal  standards  has  been  shown  to  give 
approximately  two  fold  increases  in  precision.  The  use  of 
an  argon  line  as  an  internal  standard  proved  ineffective 
and  a  measurable  inverse  correlation  was  observed  between 
analyte  and  argon  emission  signals. 

The  use  of  internal  standardisation  appears  to  warrant 
further  study.  The  construction  of  a  system  in  which 
dual  monochromators  view  the  same  segment  of  the  plasma 
(by  means  of  a  beam  splitter)  may  be  useful,  for  more 
accurately  comparing  the  performance  of  internal  standards. 

A  more  complete  investigation  of  various  other  elements 
and  emission  lines  should  reveal  which  excitation  character¬ 
istics  are  important  to  match  when  using  the  internal 
standard  principle  with  the  inductively  coupled  plasma. 
Additionally,  the  effect  of  integration  time  and  the  noise 
power  spectra  of  the  ratioed  signals  would  reveal  the 
effectiveness  of  the  internal  standard  in  removing  the 
1/f  noise  components.  Indeed,  if  the  use  of  an  internal 
standard  reduces  the  noise  to  where  shot  noise  becomes 
limiting,  then  the  use  of  long  integration  times  will  bring 
about  further  increase  in  the  precision  of  ICP  analysis. 
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APPENDIX  A 


USE  OF  ADDED  RANDOM  NOISE  TO  IMPROVE  BIT  RESOLUTION 
IN  DIGITAL  SIGNAL  AVERAGING 

Digital  signal  averaging  is  one  of  the  most  common 
instrumental  methods  in  use  today  to  improve  the  signal- 
to-noise  ratio  of  a  measurement.  It  is  important  to  keep 
in  mind  that  for  digital  signal  averaging  to  be  effective 
a  certain  amount  of  noise  must:  be  present  on  the  signal. 

In  the  absence  of  noise,  an  averaged  signal  will  be  dis¬ 
torted  by  quantization  effects  which  result  from  the 
finite  bit  resolution  of  the  analog-to-digital  converter 
(ADC)  associated  with  the  digital  signal  averager,  and  no 
amount  of  averaging  can  removethese  effects  or  improve 
measurement  precision  or  accuracy.  If  a  relatively  small 
amount  of  random  noise  is  present  on  the  signal  (standard 
deviation  of  the  noise  equal  to  about  one  half  of  the 
quantization  level,  i.e.  least  significant  bit,  of  the  ADC) 

then  the  net  effect  of  quantization  is  to  increase  the 

2 

variance  of  the  signal  by  q  /12  where  q  is  the  magnitude 
of  the  quantization  level  (A1-A3) .  Thus,  if  quantization 
effects  are  severe  (i.e.  the  original  signal  is  noise  free 
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relative  to  the  bit  resolution  of  the  ADC)  it  will  be 
advantageous  to  add  random  noise  to  the  signal  before 
averaging  the  digitized  values  (A4-A6) .  In  this  appendix 
a  simple  practical  equation  is  presented  that  can  be  used 
to  estimate  the  effective  bit  resolution  that  can  be 
achieved  by  signal  averaging  in  the  presence  of  added 
random  noise. 

As  mentioned  above,  with  sufficient  noise  on  a  signal 

the  effect  of  quantization,  under  signal  averaging,  is  to 

2 

increase  the  variance  of  the  signal  by  q  /12.  Thus  total 
signal  variance  lVT)  may  be  expressed  as: 


(qf )  2 


(Al) 


where  q  is  the  magnitude  of  the  quantization  level  and 
2 

(qf)  is  the  variance  of  the  added  noise  on  the  signal 
expressed  in  units  of  the  quantization  level.  Thus,  in 
the  situation  in  which  the  standard  deviation  (or  root-mean- 
square  value)  of  the  added  random  noise  is  twice  the 
magnitude  of  the  quantization  level,  it  would  be  2.  Signal 
averaging  in  the  presence  of  random  noise  improves  the 
signal-to-noise  ratio  [defined  as  ( signal/standard 
deviation) ]  in  proportion  to  the  square  root  of  the  number 
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_  1/2 
of  signal  repetitions  averaged  (N  ) ,  thus  the  final 

signal-to-noise  ratio  (SNR)  may  be  expressed  as: 

r 


(SNR) 


F 


(A2 ) 


where  S  is  the  magnitude  of  the  signal.  It  is  instructive 
to  express  Equation  (A2)  in  terms  of  effective  bit  res¬ 
olution  °f  the  final  signal-to-noise  ratio,  i.e. 

log2  [ (SNR) p] • : 


^BR)eff  1°g2(q)  2log2^12+f  ^  +  2log2N 


(A3) 


The  first  term  of  Equation  (A3)  represents  the  bit 
resolution  of  a  single  measurement.  It  is  important  to 
emphasize  that  this  is  not  equal  to  the  bit  resolution  of 
the  ADC  but  is  proportional  to  the  relative  magnitude  of 
the  signal  to  the  quantization  level  of  the  ADC.  The 
full  bit  resolution  of  an  ADC  is  only  utilized  when  the 
signal  fills  the  dynamic  range  of  the  ADC.  The  second 
term  represents  the  loss  in  bit  resolution  due  to  quantizing 
effects  and  the  level  of  added  noise.  The  third  term 
represents  the  build-up  of  bit  resolution  that  results 
from  signal  averaging  when  the  least  significant  bit  of 
the  ADC  is  randomized. 
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A  simple  microprocessor  controlled  data  logger  (A7) 
was  used  to  obtain  data  to  test  Equation  (A3) .  This 
system  is  described  in  detail  in  Appendix  B  and  was  based 
on  the  Intel  8748.  It  contained  an  8  bit  successive 
approximation  ADC  with  a  dynamic  range  of  0-10  V  and  a 
quantization  level  of  39  mV.  The  ADC  was  front-ended  by 
a  sample-and-hold  amplifier.  Eight  bit  conversions  could 
be  acquired  at  a  rate  of  4.4  kHz  and  up  to  65,536  suc¬ 
cessive  conversions  could  be  summed.  The  sum  was  ac¬ 
cumulated  with  triple  precision  (i.e.  24  bits).  Finally, 

32  repeat  measurements  of  the  summed  conversions  were 
automatically  acquired  and  logged  to  a  teletype  by  this 
system. 

Input  dc  level  analog  test  signals  were  generated 
using  a  16  bit  digital-to-analog  converter  (Data  Systems 
DAC-169) .  Using  this  DAC  voltage  levels  over  the  range 
0-10  V  could  be  generated  that  could  be  continuously  altered 
in  steps  as  small  as  0.153  mV.  This  DAC  was  used  to 
provide  a  stable  test  signal  to  which  various  levels  of 
random  noise  were  added.  A  General  Radio  1390B  random 
noise  generator  was  used  as  the  noise  source.  The  DAC  and 
noise  signals  were  summed  together  with  an  operational 
amplifier. 

Test  signals  were  generated  consisting  of  a  5  volt  dc 
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level  upon  which  was  superimposed  20  mV,  40  mV  and  80  mV 
rms  random  noise  (5  Hz  to  20  kHz  bandwidth) .  These  noise 
levels  correspond  to  f  =  0.5,  1  and  2.  The  data  logger 
was  used  to  acquire  32  repeat  measurements  of  one  ADC 
conversion  and  32  repeat  measurements  of  10,  100,  1000, 

10,000  and  65,536  summed  ADC  conversions.  Based  on  the  32 
repeat  measurements  the  standard  deviation  (a)  of  the 
acquired  signal  in  each  case  could  be  calculated  and  the 
signal-to-noise  ratio  (S/a)  evaluated.  The  resulting 
values  expressed  in  terms  of  bit  resolution  (log^CS/a)) 
are  presented  in  Table  AI  along  with  those  calculated 
from  Equation  (A3) .  The  agreement  is  quite  good.  Thus, 
if  the  measurement  situation  is  amenable  to  signal  averaging, 
the  effective  bit  resolution  of  an  ADC  can  be  significantly 
improved  if  averaging  is  carried  out  in  the  presence  of 
sufficient  random  noise.  Equation  (A3)  is  useful  in 
estimating  the  benefit  to  be  gained. 

This  improvement  in  bit  resolution  is  further  il¬ 
lustrated  in  Figure  Al.  The  DAC  was  again  used  to  provide 
a  test  signal,  but  this  time  the  output  of  the  DAC  was 
set  at  2.500  V  and  then  incremented  in  2.5  mV  intervals 
to  2.550  V.  The  data  logger  was  used  to  measure  1000 
summed  conversions  at  each  incremental  value,  with  and 
without  added  random  noise  (20  mV  rms  E  f  =  0.5).  Without 
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TABLE  AI 


Bit  Resolution  of  ADC  Measurements 


Conversions  Summed  Effective  Precision  of  ADC  (BITS) 

f=0 . 5  EXPERIMENTAL  EQUATION  A3 


1 

7.6 

6.8 

10 

8.5 

8.5 

100 

9.9 

10.1 

1000 

12.1 

11.8 

10000 

13.8 

13.4 

65536 

14.4 

14 . 8 

o 

• 

I — 1 

II 

4-) 

1 

6.1 

5.9 

10 

7.6 

7.6 

100 

9.1 

9.2 

1000 

11.1 

10.9 

10000 

12.5 

12.6 

65536 

13.6 

13.9 

o 

• 

CM 

II 

M-l 

1 

4.5 

5.0 

10 

6.7 

6.7 

100 

8.6 

8.3 

1000 

10.1 

10.0 

10000 

11.5 

11.6 

65536 

12.7 
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added  random  noise  the  8-bit  ADC  of  the  data  logger  can 
clearly  only  resolve  signal  changes  of  at  least  39  mV 
(i.e.  the  value  of  the  least  significant  bit)  even  if 
digitized  values  are  summed.  From  Table  AI ,  with  20  mV 
rms  added  random  noise  and  1000  summed  conversions  the 
effective  bit  resolution  should  improve  to  about  12-bits, 
making  the  least  significant  bit  equivalent  to  about  2.5  mV. 
It  can  be  seen  in  Figure  Al  that  12-bit  resolution  is 
achieved  because  2.5  mV  incremental  changes  in  DAC  output 
are  accurately  measured  when  averaging  is  carried  out  in 
the  presence  of  this  level  of  added  random  noise. 

The  examples  presented  were  based  on  an  8-bit  ADC 
because  this  is  the  most  convenient  word  length  to  man- 
iuplate  in  the  interface  to  8-bit  microprocessors.  In 
addition  some  manufacturers  have  developed  microprocessors 
with  self-contained  8-bit  ADC's  and  DAC ' s .  Examples  are 
the  Intel  8022  and  2920.  With  the  considerations  dis¬ 
cussed  here  one  need  not  think  of  these  systems  as  limited 
to  8-bit  precision  or  resolution  when  used  for  data 
acquisition.  All  considerations  presented  are,  of  course, 
equally  applicable  to  higher  resolution  conversions.  In 
fact  the  incremental  cost  in  moving  from  a  12-bit  converter 
to  a  16-bit  converter  is  quite  high.  Thus,  if  data 
acquisition  time  allows,  signal  averaging  can  be  used  to 
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achieve  the  desired  higher  bit  resolution.  Finally,  the 
considerations  presented  in  this  appendix  lead  to  some 
interesting  conclusions  about  lower  resolution  converters. 
If  the  signal  is  noisy  enough  one  should  be  able  to  relax 
the  bit  resolution  specification  of  the  ADC.  In  fact  with 
very  noisy  signals  it  is  possible  to  relax  your  specifica¬ 
tions  all  the  way  to  a  cne  bit  converter,  i.e.  a  comparator 
or  threshold  detector.  This  measurement  approach  has 
recently  been  discussed  in  the  literature  (A8)  and  such 
one  bit  converters  have  been  successfully  utilized  to 
facilitate  the  digital  data  acquisition  of  repetitive 
noisy  signals  at  very  high  sampling  rates  (A9) . 
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APPENDIX  B 


DATA  LOGGING  SYSTEM  BASED  ON  THE  INTEL  8748 
SINGLE  CHIP  MICROCOMPUTER 

One  of  the  most  common  data  logging  tasks  in  the 
laboratory  is  the  acquisition  of  repeat  measurements  of  an 
averaged  or  integrated  signal.  A  versatile  way  to  in¬ 
tegrate  a  signal  is  through  the  summation  of  a  large  number 
of  rapidly  acquired  digitized  signal  values  (Bl) .  For 
example,  if  a  measurement  system  is  capable  of  acquiring 
digitized  signal  values  at  a  rate  of  1000  conversions  per 
second,  it  is  simply  necessary  to  sum  10,  100,  1000,  or 
10,000  conversions  to  obtain  precise  signal  integration 
periods  ranging  from  10  msec  to  10  sec.  A  simple  micro¬ 
processor  based  system  was  developed  to  acquire  integrated 
values  in  this  fashion. 

The  microprocessor  chosen  for  this  application  was 
the  Intel  8748.  This  is  a  so  called  single  chip  micro¬ 
computer  in  that  it  contains,  in  one  integrated  circuit, 
all  the  main  functional  units  of  a  complete  computer  (B2) , 
i.e.  a  central  processing  unit  (CPU),  input/output  (I/O) 
ports,  a  timer/event  counter,  random  access  memory  (RAM) 
and  read  only  memory  (ROM) .  The  specific  features  of  the 
8748  are  summarized  in  Table  BI  and  a  block  diagram  of 
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TABLE  BI 


Summary  of  8748  Features 


1.  Full  8-bit  central  processing  unit 

2.  Data  memory  (RAM)  -  64  words  *  8  bits 

3.  Program  memory  (EPROM)  -  1024  x  8  bits 

4.  Input/Output  Ports 

(a)  Ports  1  and  2 

-  each  8  bits  wide 

-  all  pins  independently  programmable  as  either 
input  or  output 

-  outputs  latched 

-  inputs  not  latched 

(b)  Bus  Port 

-  8  bit  bidirectional  port 

-  all  pins  either  input  or  output 

-  outputs  latched 

-  inputs  not  latched 

5.  Input  lines 

-  3  input  lines  TO,  T1  and  testable  with  the  conditional 
jump  instruction  and  INT  an  interrupt 

6.  Timer/Counter 

-  8  bit  binary  presettable  counter 
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the  chip  is  shown  in  Figure  B1 . 

The  8748  is  one  chip  in  the  "8048"  family  of  single 
component  8-bit  microcomputers  (B3).  A  summary  of  the 
"8048"  family  of  microcomputers  is  presented  in  Table  BII. 
The  8748  is  the  development  chip  for  the  8048.  The  dif¬ 
ference  between  the  two  chips  is  that  the  8748  has  erasable 
programmable  read  only  memory  (EPROM)  while  the  8048  is 
strictly  ROM  based,  i.e.  it  must  be  factory  programmed. 

In  general,  the  8048  family  was  developed  for  high  volume 
control  applications  and  once  a  user  has  developed  a  set 
program  using  the  8748  or  the  8035  or  8039  with  off  chip 
RAM  or  EPROM,  factory  mask  programmed  8048' s  or  8049' s  can 
be  manufactured.  Within  this  family  the  8022  (see  Table 
BII)  seems  almost  ideal  for  laboratory  data  acquisition  as 
it  contains  an  on  chip  analog-to-digital  converter  (ADC) . 
However,  this  chip  is  again  ROM  based  limiting  its  use¬ 
fulness  to  quantity  users.  Thus,  the  8748  is  the  only 
single  chip  microcomputer  in  this  family  that  can  be 
effectively  utilized  by  unit  quantity  users.  However, 
as  outlined  in  Table  BI,  it  has  excellent  features  and 
capabilities  that  should  allow  the  development  of  a  very 
compact  intelligent  data  acquisition  system  with  a  minimal 
number  of  integrated  circuits. 

The  Intel  Prompt  48  unit  was  used  for  system  develop- 
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Block  diagram  8748 
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TABLE  BII 

"8048"  Family-Single  Component  Microcomputers 

8048  -  lkx 8  ROM,  64*8  RAM,  27  I/O  lines,  8-bit  timer/counter, 

full  CPU 

8748  -  EPROM  version  of  8048 
8035  -  8048  without  on-chip  ROM 

8049  -  2kx8  ROM,  128x8  RAM,  version  of  8048 
8039  -  8049  without  on-chip  ROM 

8021  -  lkx  8  ROM,  64x8  RAM,  21  I/O  lines,  8-bit  timer/counter, 

full  CPU 

8022  -  2kx 8  ROM,  64x8  RAM,  28  I/O  lines,  on-chip  8  bit 

ADC  with  sample-and-hold  and  two  multiplexed  input 
channels 


- 


ISO 


ment,  program  debugging,  and  EPROM  programming.  It  should 
be  noted  that  the  8048  family  has  its  own  language  that 
is  not  compatible  with  8080/85' s. 

A  schematic  diagram  of  the  data  acquisition  system 
is  shown  in  Figure  B2.  It  consists  of  a  sample-and-hold 
amplifier,  an  analog-to-digital  converter  and  the  8748 
microprocessor.  The  teletype  interface  for  the  8748  is 
shown  in  Figure  B3.  These  circuits  were  breadboarded 
using  a  Model  ADD-8000  analog-digital  designer  available 
from  E  &  L  Instruments.  This  unit  provides  SK-10  bread¬ 
board  cards  for  circuit  construction  which  fit  into  a 
mainframe  containing  power  supplies,  logic  level  switches 
and  indicators,  a  voltage  reference  source  and  a  function 
generator . 

The  sample-and-hold  amplifier  (AD582)  and  the  analog- 
to-digital  converter  (AD571)  were  obtained  from  Analog 
Devices.  The  sample-and-hold  amplifier  has  a  combined 
aperture  and  settling  time  (to  0.01%)  of  about  0.6  ysec 
and  an  acquisition  time  of  about  6  ysec.  The  AD571  is  a 
10  bit  successive  approximation  ADC  with  a  conversion  time 
of  25  ysec.  Only  the  eight  most  significant  bits  were 
used  in  this  application.  This  significantly  simplified 
both  the  hardware  aspects  of  data  acquisition  and  software 
complexity  which  was  particularly  important  as  programming 
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Schematic  diagram  of  data  acquisition  system. 
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space  was  at  a  premium.  An  8-bit  conversion  provides 
adequate  resolution  when  a  large  number  of  conversions  are 
summed  to  provide  an  integrated  signal  value.  This  point 
is  discussed  and  illustrated  in  Appendix  A  and  reference 
(B4) .  The  dynamic  range  of  the  ADC  was  set  at  0-10  V 
resulting  in  a  quantization  limit  of  39  mV.  The  eight 
digital  lines  from  the  ADC  were  interfaced  to  the  8748 
through  Port  1. 

The  sample-and-hold  amplifier  and  the  analog-to-digital 
converter  were  interconnected  as  recommended  in  literature 
supplied  by  Analog  Devices  (B5) .  A  conversion  is  initiated 
by  a  positive  pulse  generated  at  bit  5  of  Port  2.  This 
puts  the  sample-and-hold  amplifier  in  the  hold  mode  and 
triggers  the  ADC.  The  amplifier  is  held  in  the  hold  mode 
by  the  DR  (NOT  data  ready)  line  of  the  ADC  which  remains 
high  during  conversion.  The  return  of  the  DR  to  the  low 
state  at  the  end  of  conversion  is  sensed  at  the  TO  input 
line  of  the  8748  and  the  8-bit  value  is  read  into  the 
8748.  At  the  same  time  the  sample-and-hold  amplifier 
returns  to  the  track  mode. 

Inside  the  8748  the  new  8-bit  value  is  added  to  the 
sum  of  the  previous  conversions.  The  addition  is  carried 
out  in  triple  precision  (24  bits) .  This  allows  a  maximum 
of  65,536  full  scale  8-bit  ADC  conversions  to  be  ac- 
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cumulated.  After  the  addition,  the  program  checks  to  see 
if  the  requested  number  of  conversions  have  been  summed, 
and  if  not,  initiates  another  conversion.  The  number  of 
conversions  to  be  summed  in  order  to  achieve  the  desired 
signal  integration  period  is  hardware  set  at  pins  of 
Port  2  which  are  read  at  the  start  of  the  data  acquisition 
program.  Application  of  a  positive  1  logic  level  at  the 
appropriate  pin  selects  1,  10,  100,  1000,  10000  or  65,536 
conversions  to  be  summed.  When  the  selected  number  of 
conversions  has  been  summed  the  24-bit  binary  sum  is 
converted  to  an  8-digit  BCD  number  and  output  to  the 
teletype.  The  basic  data  rate  for  conversion  was  4.4  kHz, 
thus  a  summation  of  65,536  conversion  represents  an  in¬ 
tegration  period  of  about  15  sec.  Also  the  acquisition  of 
replicate  integrated  measurements  Was  programmed  into  the 
software.  For  our  applications  32  replicate  integrated 
measurements  were  automatically  sequenced  and  logged  to 
the  teletype.  This  enabled  the  reliable  calculation  of 
statistical  quantities  of  the  measured  values  such  as 
standard  deviation  and  signal-to-noise  ratios.  Unfortun¬ 
ately,  this  simple  8748  system  was  not  large  enough  to  do 
such  calculations  directly. 

The  8748  was  interfaced  to  the  teletype  by  means  of 
a  universal  asynchronous  receiver  transmitter  (UART )  (B6) 
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via  the  bus  port,  Port  0.  See  Figure  B3.  The  output 
routine  converts  the  BCD  digits  to  7-bit  ASCII  code  which 
the  UART  then  transmits  in  a  serial  manner  to  the  teletype. 
The  transmission  of  a  character  is  initiated  by  ap¬ 
plication  of  a  positive  pulse  to  the  data  strobe  input  of 
the  UART  from  bit  4  of  Port  2  on  the  8748.  The  UART  must 
be  externally  clocked  at  a  rate  appropriate  for  the  device 
to  which  data  is  being  sent.  Teletypes  operate  at  110  baud 
(110  bits/sec) .  The  clock  frequency,  provided  by  the 
function  generator  on  the  ADD— 8000  breadboard  system,  was 
set  at  1760  Hz.  This  is  the  correct  input  frequency  in 
order  to  allow  the  UART  to  transmit  data  at  110  baud. 

When  transmission  of  a  single  character  is  complete  the 
end  of  transmission  (EOT)  line  of  the  UART  is  sensed  at 
the  Tl  input  line  of  the  8748  and  the  next  character  is 
output.  The  rest  of  the  circuitry  in  Figure  B3  provides 
the  20  ma  current  loop  used  to  interface  the  UART  to  the 
teletype  (B6) . 

In  this  application  the  teletype  was  used  as  a  passive 
printer  and  information  could  not  be  transmitted  to  the 
8748  from  the  teletype.  One  practical  reason  was  the 
fact  that  the  unexpanded  8748  did  not  contain  sufficient 
programming  space  to  accommodate  the  necessary  programs. 

On  the  other  hand,  this  was  also  in  keeping  with  the  basic 
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philosophy  of  keeping  the  system  small.  We  did  not  want 
a  system  that  required  inputs  from  an  expensive  device 
such  as  a  teletype  in  order  to  function.  That  is  also 
why,  in  fact,  a  simple  switch  register  was  used  to  select 
the  integration  period. 

In  conclusion,  the  8748  single  chip  microcomputer  can 
be  used  to  develop  simple  intelligent  data  logging  systems. 
The  system  discussed  and  illustrated  here  consists  of  only 
four  integrated  circuits,  the  sample-and-hold  amplifier, 
the  analog-to-digital  converter,  the  8748  and  the  UART .  It 
is  capable  of  acquiring  digitized  signal  values  with  a 
resolution  of  8-bits  at  rates  up  to  4.4  kHz.  Pre-selected 
numbers  of  conversion  can  be  summed  to  obtain  integrated 
signal  measurements.  Summation  is  carried  out  in  triple 
precision  (24  bits)  and  up  to  32  replicate  measurements 
can  be  sequentially  acquired  and  logged,  in  BCD  format,  to 
a  serial  terminal  such  as  a  teletype. 
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APPENDIX  C 


SOFTWARE  DEVELOPED  FOR  DATA  ACQUISITION  AND  PROCESSING 

The  software  for  data  acquisition  and  processing  listed 
here  are  written  in  Fortran  IV  and  Macro  assembly  language; 
for  implementation  on  PDP  11/10  minicomputers.  While  the 
majority  of  the  programs  can  be  run  on  computers  with 
16  k  bytes  core,  the  Fourier  transform  routine  requires  a 
32  k  bytes  core. 

The  programs  are  documented.  Program  RMBINT  allows 
background  subtracted  signals  to  be  acquired  from  a  photo¬ 
multiplier  tube  and  integrated  for  specified  periods.  This 
program  uses  sub-routine  RBPROG  written  in  assembly 
language  to  actually  acquire  and  integrate  the  signal. 

The  program  RMBPMT  is  part  of  the  noise  power  spectra 
package.  It  allows  time  series  data  to  be  collected  from 
a  photomultiplier  tube.  A  specified  number  of  scans  may 
be  collected  so  that  signal  averaging  may  be  carried  out 
after  the  Fourier  transformation  step.  Sub-routine  AQUIRE 
written  in  assembly  language  is  used  for  the  actual  data 
acquisition  step.  This  sub-routine  was  written  for 
simultaneous  dual  channel  data  acquisition  using  external 
sample-and-hold  amplifiers.  However,  in  this  particular 
application  only  one  channel  is  utilised. 
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After  acquisition  by  program  RMBPMT  data  is  d.c. 
level  subtracted,  apodised  and  Fourier  transformed  by 
program  RBFFTP .  This  program  also  provides  signal 
averaging  ability,  after  the  Fourier  transformation. 

Noise  power  spectra  are  finally  plotted  using  the 
plotting  routine  RBPLOT . 

Dual  channel  data  acquisition  for  correlation  studies 
was  achieved  by  the  use  of  program  CORPMT  and  sub-routine 
AQUIRE.  CORPMT  not  only  allows  dual  channel  data  to  be 
acquired  simultaneously  but  also  calculates  the  correlation 
coefficient  between  the  two  sets  of  data. 
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C  RMBINT 

C  PMT  SIGNAL  INTEGRATION  ROUTINE 

C 

C  FILES  REQUIRED 

C  1 >  RBPROG 

C  2)  131  BLOCK  FORTRAN  LIBRARY  (FORLIN) 

C 

C 

DIMENSION  SIGN( 100) t BSB ( 1 00 ) * STDEV ( 100 ) , SNR ( 1 00 ) 

DOUBLE  PRECISION  ARRAY ( 100 >» RRPTS . SUM • SMSQ 

COMMON  /RON/JINT * IBL» IBH 

MNEX=0 

1000  MNEX=MNEX+1 

WRITE ( 7 » 1 4  ) 

14  FORMAT (  '  'INPUT  INTEGRATION  TIME.') 

READ (5» 16) TIME 
16  FORMAT ( F6 . 2 ) 

WRITE (7*18) 

18  FORMAT < '  '  *  '  DO  YOU  REQUIRE  BACKGRND  SUBTRACTION  Y=1  N=0  .') 

READ(5*20) IBACK 
20  FORMAT  (ID 

WRITE(7f 22) 

22  FORMAT ( '  '» 'INPUT  NUMBER  OF  REPEATS.') 

READ(5*24)MRPTS 
24  FORMAT ( 13 ) 

IF( IBACK. EQ.O)  GO  TO  1006 

1005  WRITE (7  *26 ) 

26  FORMAT ('  ' r 'ENTER  0  FOR  BACKGROUND  1  FOR  SAMPLE.') 

READ (5, 20) I  SOB 

1006  JINT=IFIX(TIME*1000.0) 

NRPTS=MRPT5 

C 

C  CALLS  MACRO  DATA  ACQUISITION  ROUTINE 

C 

WRITE ( 7 » 28 ) 

28  FORMAT ( '  '» 'ENTER  CR  TO  START  AQUISITION  .  '  ) 

READ (5 f 20) ISTRT 
1015  CALL  RBCONU 

C 

C  CONVERTS  TWO  WORD  INTEGER  DATA  TO  FLOATING  POINT 

C 

ARRAY < NRPTS )=DBLE( FLOAT ( I BL ) >+DBLE( FLOAT (IBH) ) *32768.0 
C 

C  OUTPUTS  DATA  AS  IT  ARRIVES 

C 

N=MRPTS— NRPTS+1 

WRITE ( 7 » 66 ) N  >  ARRAY (NRPTS ) 

66  FORMAT (I5f5XfE16.7) 

NRPTS=NRPTS-1 
IF ( NRPTS . NE . 0 )  GO  TO  1015 
C 

C  THIS  SECTION  CALCULATES  MEAN  AND  STD. DEV.  OF  ARRAY 

C 

SUM=0 . 0 

DO  100  NRPTS=1,MRPTS 
SUM=SUM+ARRAY( NRPTS) 

100  CONTINUE 

SMSQ=0 . 0 


DO  200  NRPTS=1 rMRPTS 
SMSQ=SMSQ+ ARRAY  <  NRPTS  >  #*2 
200  CONTINUE 

RRPTS  =  DBLE( FLOAT ( MRPTS  > ) 

AURG=SNGL ( SUM/RRPTS ) 

STDEU ( MNEX )  =SNGL  ( DSCJRT  (  (  SMSCJ-SUM**2/RRPTS  )  /  ( RRPTS-1 . 0  )  )  ) 

C 

C  THIS  SECTION  SORTS  OUT  SIGNAL  OR  BACKGROUND  DATA 

C 

IF ( IBACK i EQ . 0 )  GO  TO  1020 
IF ( I SOB . EQ . 1 )  GO  TO  1025 
UR I TE  <  7  >  30 ) AVRG  » STDEU ( MNEX ) 

30  FORMAT ( '  '»10H  BACKGRND= , E10 . 4 , 1 2H  STD-DEU= • E10 , 4 > 

BANVAL=AVRG 
GO  TO  1005 
1020  BANUAL=0 . 0 

1025  SI GN ( MNEX ) =AURG 

BSB(MNEX)=SIGN(MNEX)-BAKUAL 
SNR  (  MNEX  )  =BSB  ( MNEX  >  /STDEU  ( MNEX  > 

URITE  <  7f  32  >  MNEX  *  BSP  (  MNEX  >  *  STDEV MNEX )  » SNR  <  MNEX  ) 

32  FORMAT <  '  ' r 13 » IX » ' S IGNAL= ' , E 1 0 , 4 , 4X f 'STD-DEU=' fEIO ,4f4Xf 'SNR='  • 

1  E10.4) 

1027  URITE  <7»  34 ) 

34  FORMAT (  '  'f'TYPE  1  TO  EXIT  ,2  TO  RUN  PROGRAM  AGAIN  -3  TO  LIST ' > 

READ(5f20) I RUN 
GO  TO  <1040f1000f1 300 )  IRUN 
1030  URITE <  7  f  36 ) 

36  FORMAT  (  '  '  f'  SIGNAL  BSSIGNAL  STD-DEU  SNR'/) 

DO  700  NEX=1 1  MNEX 

URITE(7f38)NEXfSIGN<NEX) f BSB ( NEX ) f  STDEU ( NEX ) f SNR < NEX) 

38  FORMAT ( '  ' » 13 f 3X f E 1 0 , 4 , 3X f E10 . 4 , 3X . E? . 3 r 3Y . E? , 3 ) 

700  CONTINUE 

GO  TO  1027 
1040  STOP 

END 


=s 


.TITLE  RBPROG 
. MCALL  » • V2 . «  .  . REGDEF 
. .V2. . 

.REGDEF 
. GLOBL  RBCONV 

ADSR=170400 
ADBF=1 70402 
CKSR=170404 
CNPR=170406 


RBCONV: 

MOV 

JINT  » RO 

9 

MOV 

*176027  . CKPR 

MOV 

*000403. CKSR 

MOV 

♦000000 • R2 

MOV 

♦000000 . R3 

LP2 : 

MOV 

*000040. ADSR 

L.  pi  : 

TSTB 

ADSR 

BPL 

LP1 

MOV 

ADBF.R4 

ADD 

R4.R2 

ADC 

R3 

DEC 

RO 

BNE 

LP2 

9 

ASL 

R3 

ROL 

R2 

ADC 

CLC 

R3 

ROR 

R2 

MOV 

R2> IBL 

MOV 

R3. IBM 

( SET  CONVERSION  COUNTER 

(RATE  TO  CLK  PRESET  REGISTER 
» SET  CLK  STSTIJS  REGISTER 
( CLEAR  REGISTER  2 
( CLEAR  REGISTER  3 
(START  CONVERION  CHO 
( CHECK  FOR  EOC 
(  LOOP 

( MOVE  DATA  TO  REGISTER  4 
(  PERFORM  INTEGRATION 
(TAKE  CARE  OF  OVER  FLOW 
(DECREMENT  CONVERSION  COUNTER 
(LOOP  FOR  NEXT  CONVERSION  IF  REOD 

(CONVERTS  DATA  TO  FORTRAN  I  FORMAT 


jint: 

ibl: 

IBHJ 


.EVEN 

. CSECT  RON 
.WORD  0 
.WORD  0 
.WORD  0 
.CSECT 
.EVEN 
.END 
* 


END 


RBCONV 


' 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

20 

5010 

5020 

5030 

5040 


5050 


50<f>0 


5070 

C 

C 

c 

c 

c 

c 


c 

c 

c 


5080 


c 

c 

C 


RMBPMT  •  FOR' 

DATA  ACQUISITION 

SUBROUTINES  NEEDED J 

( 1 )  AQUIRE 

( 2 )  SYSL IB 

( 3 )  FORLIN 


DIMENSION  TITLE ( 3 )  f  DAT ( 3 > 

COMMON  /EON/  JNP  >  JSR  t  Ji  *R  •  I A  ( 1032  )  /RMB/  IB  (1032) 

INPUT  ACQUISITION  PARAMETERS 

CALL  DATE ( DAT ) 

URI TE ( 7 » 5010 ) 

FORMAT ( '  ENTER  NUMBER  OF  DATA  SETS  REQUIRED') 

READ ( 5  1 5020 )  NFILE 
FORMAT (14) 

UR  I TE ( 7 » 5030 ) 

FORMAT ( '  '>'  INPUT  SAMPLE  NAME') 

READ ( 5 1 5040 )  TITLE 
FORMAT (3A4) 

URITE ( 7 1 5050 ) 

FORMAT ( '  '•'  INPUT  NUMBER  OF  POINTS') 

READ ( 5  f 5020  )  JNP 
URITE ( 7  > 5060 ) 

FORMAT ('  '» 'INPUT  CLOCK  RATE ' t / r '  2=1  MHZ'r/r'  4=100  KHZ'./f 
6=10  KHZ ' » / » '  8=1  KHZ'f/r'  10=100  HZ'f/f'  14=60  HZ'/) 

READ ( 5  f  5020 )  IMCR 
WRITE ( 7  t 5070 ) 

FORMAT ( '  INPUT  NUMBER  OF  CLOCK  PULSES  BETWEEN  SAMPLES') 

READ  <  5  f  5020 )  JSR 

SET  UP  LPS  CLOCK  STATUS  REGISTER  PARAMETERS 
JCR=( IMCR+257) 

CALCULATE  ACQUISITION  FREQUENCY 

IF  (IMCR.EQ.2)  FREQ= 1000000 . O/FLOAT  <  JSR ) 

IF  (IMCR.EQ.4)  FREQ=  1 00000 . 0/FL.0AT  (  JSR  ) 

IF  (IMCR.EQ.6)  FREQ=10000 . O/FLOAT ( JSR  ) 

IF  (IMCR.EQ.8)  FREQ=1000 . O/FLOAT ( JSR ) 

IF  (IMCR. EQ. 10)  FREQ= 100 . O/FLOAT ( JSR ) 

IF  (IMCR. EQ. 14)  FREQ= 60 . O/FLOAT ( JSR ) 

IFREQ=IFIX ( FREQ ) 

ACQUIRE  DATA  UITH  SUBROUTINE  AQUIRE 
WRITE ( 7  f  5080 ) 

FORMAT ( '  ','  CR  TO  START  ACQUISITION') 

READ ( 5 1 5020 )  I 
DO  500  IFILE*1 »  NFILE 
CALL  AQUIRE 

FORMS  TEMPORARY  DATA  FILE  ON  DISC 


. 


n  n  n 


JFILE=IFILE+20 
I A ( JNP+ 1 ) = JNP 
I A ( JNP+2 ) =IFREQ 
I A  <  JNP +3  >  =  1 A ( 1 ) 

I A ( JNP  +  4 ) = 1 
I A ( JNP+5 )  =  1 A  ( 1 ) 

I A ( JNP+6 ) =1 

CALL  ASSIGN  <  JFILE ) 

WRITE  <  JFILE  >  < IA< I) »I*1 » JNP) 

ENDFILE  JFILE 
CALL  CLOSE ( JFILE ) 

500  CONTINUE 

C 

C  PROGRAM  CONTROL 

C 

90  WRITE (7» 5090 ) 

50°0  FORMAT ( '  '»///'  STORE  DATA  =1*  RERUN  =2f  EXIT  =3') 

READ  <  5 » 5020 )  J 
GO  TO  ( 10 » 20  r 30 )  J 

COMPOUNDS  DATA  FILES  *  FORMS  PERMANENT  DATA  FILE  ON  DISC 
10  WRITE ( 7  j 4000 ) 

6000  FORMAT ( '  ENTER  DEI  FILE . DAT  FOR  STORAGE'//) 

CALL  ASS  I GN ( 2 • ' DK t FTN2 . DAT ' , -1 f 'NEW' » ' NC  '  • 1 ) 

WRITE ( 2 )  NFILE 
WRITE ( 2 )  JNP 
WRITE ( 2  )  IFREQ 
WRITE (2 ) ( DAT (I)»I=1»3) 

WRITEC2) ( TITLE ( I >  » 1  =  1 r3> 

DO  600  JFILE=21 -NFILE+20 
CALL  ASSIGN( JFILE) 

READ (JFILE) ( IA ( I ) *  1  =  1 > JNP) 

WRITE (2) ( IA( I ) » 1=1 » JNP) 

ENDFILE  JFILE 
CALL  CLOSE ( JFILE ) 

600  CONTINUE 

ENDFILE  2 
CALL  CLOSE ( 2 ) 

GO  TO  90 
30  STOP 

END 

* 


, 


♦TITLE  TWO  CHANNEL  ADC 
♦ MCALL  ♦ ♦ V2 ♦ ♦ r ♦ REGDEF 
♦ .V2. . 

♦REGDEF 
♦ GLOBL  ARUIRE 
i 
9 

ADSR=170400 
ADBF=170402 
CKSR=170404 
CKPR=170406 
DBUF  = 1704 1 4 
DBSR=170410 


ARUIRE : 

MOV 

JNP ♦ RO 

i NUMBER  OF  POINTS  TO  RO 

CLR 

DBSR 

f CLEAR  DIGITAL  BUFFER  STATl. 

MOU 

*177777 ,DBUF 

? CAUSE  ANALOG  DATA  TRACK 

NEG 

JSR 

J TWO'S  COMPLEMENT 

MOV 

JSRfCKPR 

t RATE  TO  CLOCK  PRESET 

MOV 

JCRfCKSR 

i SET  CLOCK  STATUS  REGISTER 

MOV 

♦  I A  ♦  R  2 

i SET  ADDRESS  ARRAY  I A 

• 

9 

MOV 

♦ IB ♦ R3 

> SET  ADDRESS  ARRAY  IB 

* 

9 

LP3  ♦ 

TSTB 

CKSR 

*  CHECK  CLOCK 

BPL 

LP3 

i  LOOP 

» 

9 

{ 

i 

MOV 

JCR » CKSR 

i RESET  MODE  FLAG 

MOV 

*177776 » DBUF 

.CAUSE  ANALOG  DATA  HOLD 

CLR 

ADBF 

.CLEAR  A/D  BUFFER 

MOV 

*000001 fADSR 

» SELECT  CHANNEL  0 

lpi: 

TSTB 

ADSR 

( CHECK  FOR  ADC 

BPL 

LPI 

iLOOP 

MOV 

ADBF. (R2>+ 

f STORE  DATA  ARRAY  IA 

i 

CLR 

ADBF 

{CLEAR  A/D  BUFFER 

MOV 

*000401 .ADSR 

.SELECT  CHANNEL  1 

LP2: 

TSTB 

ADSR 

{CHECK  FOR  ADC 

BPL 

LP2 

{LOOP 

MOV 

ADBF. (R3)+ 

{STORE  DATA  ARRAY  IB 

i 

i 

i 

MOV 

*177777. DBUF 

{CAUSE  ANALOG  DATA  TRACK 

DEC 

RO 

{DECREMENT  POINT  COUNTER 

BNE 

LP3 

{LOOP  IF  NOT  ZERO 

i 

* 

NEG 

JSR 

I  TWO ' S  COMPLEMENT 
• 

9 

* 

9 

♦  EVEN 

♦  CSECT 

RON 

jnp: 

♦  WORD 

0 

jsr: 

♦  WORD 

0 

JCR  ♦ 

♦  WORD 

0  . 
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ia: 

•  BLKU 

1032 

• CSECT 

RMB 

if: 

» BLKW 

1032 

♦  r.SECT 

.  EVEN 

end: 

.END 

AQUIRE 

* 
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c 

c 

c 

c 

c 

c 

c 

c 


1200 

1000 


1100 


2000 

2100 


1300 

1400 

1A00 

1700 

1800 

* 


c 

c 

c 


100 


RBFFTP , FOR 

FLOATING  POINT  FAST  FOURIER  TRANSFORM 

files  required: 

<1 )SYSLIB 
( 2  >  FORLI B 
( 3 )LPSL IB 

DIMENSION  X'llOO) ,Y(1100) ,L(20) 

DIMENSION  DAT'D  ,TITLE(3> , AU(20) 

INTEGER  ID (1100) 

CALL  DATE (DAT) 

UR  ITE  ( 7,1 200  )  (DAT'D  ,1  =  1,3) 

FORMAT ( '  ' , 3A4 ) 

WRITE (7, 1000 ) 

FORMAT ('  ', 'ENTER  DE? FILE: DAT  FOR  STORAGE'/'  ' ) 

CALL  ASSIGN ( 12, 'DKJFTN12.DAT' ,-l , ' NEW ' , 'NC' , 1 > 

WRITE (7,1 100 ) 

FORMAT ( '  ', 'ENTER  DE  *  FILE  »  DAT  FOR  DATA'/'  ') 

CALL  ASSIGN (10, 'DKtFTN10.DAT' ,-l , 'OLD' , 'NC' *1 ) 

WRITE ( 7 , 2000 ) 

FORMAT ( '  ' , ' AP0DI2ATI0N?  0  FOR  NO,  OTHERWISE  FACTOR'1 
READ ( 5 ,2100  )  IAPOD 
FORMAT (II) 

READ ( 10 )  NF 
READ ( 10 )  NP 
READ ( 10  )  NS 
READ(10)(DAT(I) ,1=1,3) 

READ' 10) (TITLE' I ) , 1=1 ,3) 

WRITE (12)  NF 
WRITE ( 12)  NP 
WRITE( 12)  NS 
WRITE' 12) ' DAT ( I ) , 1=1 ,3) 

WRITE( 12) ( TITLE ( I ) , 1=1 ,3) 

WRITE  (7,  .1.300)  ( DAT  (  T  )  ,  1=1,3) 

FORMAT''  ','DATEJ  ',3A4> 

WRITE (7, 1400) 'TITLE' I ) , 1=1 ,3) 

FORMAT (  '  ', 'TITLE:  '*3A4> 

WRITE (7,1 600 )  NF 

FORMAT (  '  ', 'NUMBER  OF  SETS  OF  DATA :  ',14) 

WRITE (7, 1700)  NP 

FORMAT (  '  ', 'NUMBER  OF  DATA  PER  SET?  ',14) 

WRITE ' 7 , 1800 )  NS 

FORMAT ( '  ', 'SAMPLING  FREQUENCY  OF  DATA  ACQUISITION  ?  ' 

'  HZ'/'  ') 

CALL  ASSIGN ( 14, 'DKtFTN14.DAT' ,12, 'SCR' , 'NC'  ,1 > 

DO  1500  IF=1 » NF 
CALL  LED( IF, ' 15' ) 

DC  LEUEL  SUBTRACTION 

READ  <10)(ID(I),I  =  1,NP) 

SUM=0 . 0 

DO  100  1=1, NP 

X ( I ) “FLOAT ( ID ( I ) ) 

SUM=SIJM  +  X  ( I ) 

CONTINUE 

SUM=SUM/FLOAT ( NP ) 


,17, 


- 
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DO  200  1  =  1  *  NP 
X  <  I ) =X ( I ) -SUM 
200  CONTINUE 

AV ( IF ) *SUM 
C 

c  APODISATION 

C 

IF ( I APOD . EQ . 0 )  GO  TO  2600 
FACTOR»FLOAT< IAPOD) 

ZPD=512,0 
DO  2500  1=1 ,NP 

DDUMO=ABS  <  FLOAT ( I ) -ZPD )/ f  FLOAT ( NP ) -ZPD ) 
DDUM1=EXP<-FACT0R*DDUM0**2) 

X(  I  >=X<  I  XDDUN1 
2500  CONTINUE 

C 

C  DC  LEVEL  SUBTRACTION 

C 

SUM=0 . 0 
DO  120  1=1. NP 
SUM=SUM+X( I ) 

120  CONTINUE 

SUH=SUM/FLOAT(NP> 

DO  220  1=1, NP 
X( I )=X( I ) —SUM 
220  CONTINUE 

C 

C  FFT  ALGORITHM 

C 

2600  NPT  =  NP/2 

N=NP 

DO  2  1=1 .NPT 
X( I ) =X( 1*2-1 ) 

2  Y(I)=X(I*2) 

DO  3  I=NPT+ 1 » NP 
Y  < I ) =0 . 0 

3  X(I)=0.0 
N2P0U=10 
NTHP0W=2**N2P0W 
N4P0U=N2F’0W/2 

IF  ( N4F'0W  )60 .63  »  60 
60  DO  61  IPASS=1 .N4P0W 

NXTLTH=2**(N2P0W-2*IPASS> 

LENGTH=4*NXTLTH 

SCAL  =  6. 283 1853/FLO AT (LENGTH) 

DO  61  J=1 , NXTLTH 
ARG=FLOAT( J-l >*SCAL 
C1=C0S ( ARG ) 

S1=S IN ( ARG ) 

C2=C1*C1-S1*S1 

S2=C1*S1+C1*S1 

C3=C1*C2-S.1*S2 

S3=C2*S1+S2*C1 

DO  61  I  SQL..OC=LENGTH  • N THROW ■ LENGTH 

J1=I SQLOC-LENGTH+ J 

J2= Jl+NXTLTH 

J3= J2+NXTL.TH 

J4= J3+NXTLTH 

R1=X ( J1 )+X( J3) 

R2=X ( J1 ) -X  <  J3 ) 

R3=X ( J2 ) +X ( J4 ) 

R4=X ( J2 ) -X ( J4 ) 

F I 1=Y ( J1 ) +Y ( J3 ) 

FI2=Y(J1)-Y(J3) 

FI3=Y(J2>+Y(J4) 

FI4=Y( J2)-Y( J4) 
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X( JI >=R1+R3 
Y ( J1 ) =FI 1+FI3 
IF ( J— 1 )64?62>64 

64  X< J3)=C1*<R2+FI4)+S1*'FI2-R4) 

Y  (.  J3  )=-Sl*(R2+FI4)+Cl*<FI2-R4) 
X< J2)=C2*<R1-R3)+S2*(FI1-FI3> 

Y ( J2 )=-S2* ( R1—R3  >  +C2* ( FI 1— FI3  > 
X< J4)=C3*(R2-FI4)+S3*(R4+FI2) 
Y< J4)=-S3*<R2-FI4)+C3¥(R4+FI2> 
GO  TO  61 

62  X(J3)=R2+FI4 
Y< J3)=FI2-R4 
X< J2)=R1-R3 

Y (J2)=FI 1— FI3 
X< J4)=R2-FI4 
Y ( J4 ) =R4+FI2 
61  CONTINUE 

63  IF ( N2P0W-2*N4P0W ) 65 » 66  » 65 

65  DO  67  J=1 » NTHF’OW  f  2 
R1=X( J)+X( J+l ) 

R2=X ( J ) -X ( J+l ) 

FI1=Y( J)+Y( J+l ) 

FI2=Y( J)-Y< J+l ) 

X< J)=R1 
Y ( J ) =FI 1 
X( J+l )=R2 

67  Y ( J+l ) =FI2 

66  DO  68  J=1 » 13 
L ( J ) =1 

I F  (  J— N2P0W  )6CJ»69*68 
6?  L< J)=2**<N2P0U+1-J) 

68  CONTINUE 
I  J=1 

L1=L ( 13 ) 

NTHPOW=L< 12) 

I SQLOC=L (11) 

J=L( 10) 

N2P1=L ( 9 ) 

N2=L ( 8 ) 

NP2M J=L ( 7 ) 

L8--L  ( 6 ) 

N2P0U=L(5) 

N4F'0W=L  (4) 

LENGTH=L ( 3 ) 

NXTLTH=L ( 2 ) 

IPASS=L< 1 ) 


DO 

601 

J 1  =  1 » L 1 

DO 

601 

J2= JI  •  NTHPOKI  ?  L 1 

DO 

601 

J3=J2 » ISQLOC  • NTHPOW 

DO 

601 

J4=J3* Jf ISOLOC 

DO 

601 

J5=J4  r N2P1 t J 

DO 

601 

J6= J5 » N2  f N2P1 

DO 

601 

J7=J6fNP2MJ»N2 

DO 

601 

J8  =  J7»L8»  NP2H J 

DO 

601 

J9  =  J8 » N2P0U  »L8 

DO 

601 

J10=J9»N4P0U>N2P0I.I 

DO 

601 

J11=J10»LENGTH*N4P0W 

DO 

601 

J12* JI 1 f  NXTLTH  f LENGTH 

DO 

601 

JI=J12t IPASSfNXTLTH 

IF< I J-JI )610»610»601 
610  R=X ( I J ) 

X(  I J ) =X ( JI ) 

X< JI )=R 
FI  =  Y  < I J ) 

Y< I J)=Y<  JI ) 

Y( JI ) =FI 


601  IJ=IJ+1 

ARG=3 . 1 415927/FLOAT ( N ) 

C1=C0S ( ARG ) 

S1=-SIN(ARG> 

C1JX=1 . 

S1JX=0. 

N2=N/2 

N2P1=N2+1 

DO  70  J=2»N2F'l 

NF'2M  J=N+2- J 

SORR 1  =X  (  J  )  +X  ( NP2M  J ) 

SORI  1  =  Y  (  J )  -  Y  (  NP2M  J ) 

R=C1 JX 

Cl JX=C1 JX*C1-S1 JX*S1 
S1JX=R*S1+S1JX*C1 
S0RR2  =  X ( J ) -X ( NP2M  J ) 

S0RI2  =  Y<  J)+Y<NF'2MJ) 

S0RR3=C1 JX#S0RR2-S1 JX*S0RI2 
SORI  3=C1  JX*S0RI2  +  S1  JY  YS0PF:2 
Y  <  J )=0 . 5*  <  SORI 1-S0RR3 ) 

X< J)=0.5*(S0RR1+S0RI3> 

IF ( J-N2P1 ) 71 ,70,71 
71  Y ( NP2M J )=0 ♦ 5# (— SORI 1-50RR3 ) 

X<NP2MJ  >=0 .  5#  < S0RR1-S0PI3 ) 

70  CONTINUE 

X( 1 )=X( 1 )+Y( 1 ) 

Y ( 1 ) =0 . 

DO  501  1  =  1 , 1024 

X(I )  =  (X( I )*X(I )+Y( I )*Y(I >  ) 

501  CONTINUE 

WRITE< 14) (X( I ) r 1=1 ,NP ) 

1500  CONTINUE 

C 

C  DISC  CONTROL 

C 

ENDFILE  10 
CALL  CLOSE (10) 

ENDFILE  14 
REWIND  14 
C 

C  SIGNAL  AVERAGING  FOURIER  DOMAIN 

C 

DO  3000  1=1, NP 
3000  Y ( I ) =0 . 0 

DO  3100  1=1, NF 
READ  <14)<X(J) , J=1 ,NP) 

DO  3200  J=1,NF 
3200  Y(J)=Y(J)+X(J) 

3100  CONTINUE 

DO  4000  1=1, NP 
Y< I >=Y< I >/NF 

4000  CONTINUE 
YAV=0 . 0 

DO  4001  1=1, NF 
YAV=YAV+AV< I ) 

4001  CONTINUE 
YAV=YAV/FLOAT ( NF ) 

ENDFILE  14 

CALL  CLOSE ( 14  ) 

C 

C  WRITE  FT  DATA  ON  DISC 

C 

WRITE < 12 >  <  Y ( I > » 1  =  1 ,NP ) 

WRITE( 12) ( YAU) 

ENDFILE  12 
CALL  CLOSE ( 1 2  > 

STOP 

END 
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c 

c 

c 

c 

c 

c 

c 

c 


1200 

1100 

3500 

6100 

6000 


1300 

1400 

1600 

1700 

1800 

*' 

1900 

C 

C 

C 

4000 


400 


RBPLOT . FOR 

ZETA  PLOTTING  FOR  FOURIER  TRANSFORMS 

FILES  required: 

(1 )SYSLIB 
< 2 ) XYL IB 
<3)F0RLIB 

DIMENSION  X(1100) *Y< 1100) 

DIMENSION  DAT ( 3 ) *  TI TLE  (  3  > 

CALL  DATE (DAT) 

WRITE (7. 1200) ( DAT ( I ) * 1=1 >3) 

FORMAT ( '  '  *  3A4  ) 

WRITE (7*1 100 ) 

FORMAT ('  '♦ 'ENTER  DE  t  FILE . DAT  FOR  DATA'/'  ') 

CALL  ASSIGN (10* ' DK \ FTN10 > DAT ' *-l * ' OLD ' *  'NC' -1) 

WRITE ( 7  *  3500 ) 

FORMAT ( '  ', 'ENTER  LENGTH  FOR  X-AXIS*  INCH') 

READ ( 5  *  6000  > I IXX 
WR ITE ( 7  *  6100 ) 

FORMAT ( '  '* 'ENTER  LENGTH  FOR  Y-AXIS*  INCH') 

READ(5*6000)IIYY 

FORMAT( 12) 

XAXIS=FLOAT (IIXX) 

YAXIS=FLOAT( IIYY) 

READ ( 10 )  NF 
READ (10)  NP 
READ ( 10  )  NS 
READ( 10) (DAT( I ) *1=1 *3) 

REAIK  10) (TITLE(I) *1  =  1*3) 

READ( 10) (Y( I ) *  1  =  1 » NP ) 

READ (10)  YAU 

WRITE ( 7  *  1300 ) ( DAT ( I ) *  1  =  1*3) 

FORMAT (  '  '*'DATEJ  '*3A4) 

WRITE ( 7 *1400) (TITLE(I) *1=1*3) 

FORMAT (  '  ' » 'TITLE  *  '*3A4) 

WRITE (7*1 600 )  NF 

FORMAT < '  ' * ' NUMBER  OF  SETS  OF  DATA:  '*14) 

WRITE(7* 1700)  NP 

FORMAT (  '  '» 'NUMBER  OF  DATA  PER  SET!  '*14) 

WRITE (7*1 800  )  NS 

FORMAT ( '  '* 'SAMPLING  FREQUENCY  OF  DATA  ACQUISITION:  '*17, 
HZ'/'  ') 

WRITE(7* 1900)  YAU 

FORMAT (  '  ' * ' MEAN  SIGNAL=  '*E11.4) 

ENDFILE  10 
CALL  CLOSE ( 10 ) 

NORMALISES  SPECTRUM  IF- REQUIRED 
WRITE (7 , 4000 ) 

FORMAT ( '  '* 'NORMALISED  SPECTRUM  =0*  NON-NORMAL ISED  =1') 

READ(5*9001)IN0RM 

IF  (INORM.EQ.l)  GO  TO  700 

DO  400  1=1 »NP 

Y ( I )=Y< I )/YAV 

CONTINUE 


c 

C  CALCULATES  S  DISPLAYS  MAX  AMPLITUDE 

C 

200  DO  500  1=1. NP 

IF  <Y< I ) .GT.AMAX)  AMAX=Y(I) 

500  CONTINUE 

UR I TE ( 7  »  8000 ) AMAX 

8000  FORMAT ( '  '.'MAX  AMPLITUDE  =  '.E10,4> 

C 

C  COMPUTES  SCALING  FACTOR 

C 

WRITE ( 7 » 6050 ) 

6050  FORMAT ( '  '.'DO  YOU  REQUIRE  SCALED  PLOT  Y=1  N=0'> 

READ ( 5 .8100  )  M 
8100  FORMAT (II) 

IF(M.EQ.O)  GO  TO  5 
Y ( NP  +  1  )=0.0 
Y (NP+2)=AMAX/YAXIS 
GO  TO  6 

5  WRITE (7.6150) 

6150  FORMAT ( '  '.'INPUT  MAX  AMPLITUDE  TO  BE  PLOTTED') 

REA D (5. 6200) P MAX 
6200  FORMAT (E10 . 4 ) 

Y ( NP+1 ) =0 . 0 
Y(NP+2)=PMAX/YAXIS 

6  RSCAL=FL0AT(NS)/2.0 
RINC=PSCAL/FLOAT(NP) 

DO  6500  1=1. NP 
J=I-1 

6500  X< I )=FLOAT< J)*RINC 

C 

C  TRUNCATES  DATA  TO  FIT  PLOT 

C 

DO  600  1=1. NP 

IF ( Y ( I ) . GT . PMAX )  Y ( I ) =PMAX 
600  CONTINUE 

C 

C  XY  PLOTTING  ROUTINE 

C 

CALL  PLOTST ( 0 . 005 » ' IN  '  ) 

CALL  PLOT  (0.2. 0,5. -3) 

CALL  SYMBOL <0,0. 0.0.0. 16. DAT. 90 .0. 12) 

CALL  SYMBOL (0.0. 2. 0.0. 16. TITLE* 90.0. 12) 

CALL  PLOT < 0 . 8 » 0 . 0  r -3 ) 

CALL  AXIS(0. 0.0.0. 5HP0WER ,+5.YAXIS.90*0.Y(NP+l ) .Y(NP+2) ) 
CALL  SCALE (X.XAXIS.NP.l ) 

CALL  AXIS(0.0»0,0» 9HFREQUENCY »— 9.XAXIS.0*0.X(NP+1 ) .X(NP+2> ) 
CALL  LINE(X.Y.NP. 1.0.0) 

CALL  PLOT (XAXIS.0.0.-3) 

CALL  PLOT (0.5. -0 ,5.-3) 

CALL  PLOTND 

STOP 

END 
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CORPMT .FOR 

TWO  CHANNEL  DATA  ACQUISITION  FOR  CORRELATION 

SUBROUTINES  NEEDED: 

( 1 )  AQUIRE 

( 2 )  SYSLIB 
( 3  >FORLIN 


DIMENSION  TITLE ( 3 ) » DAT (3) f A( 1200) »B( 1200) 

DOUBLF  PRECISION  SUM A  » SUMB  t SUMA2 r SUMB2  r SUMAB  »  SA  t SB • SAB t R 
COMMON  /RON/  JNP r JSR t JCR r I A ( 1024 )  /RMB/  IB'1024> 

INPUT  ACQUISITION  PARAMETERS 


CALL  DATE ( DAT ) 

FORMAT (14) 

WRITE  (.7 1 5030 ) 

FORMAT ('  INPUT  SAMPLE  NAME') 

READ ( 5  f 5040 )  TITLE 
FORMAT ( 3A4 ) 

WRITE(7f 5050) 

FORMAT ('  ' » '  INPUT  NUMBER  OF  POINTS') 

READ ( 5 1 5020 )  JNP 
WR I TE ( 7  ? 5060 ) 

FORMAT ( '  ' » ' INPUT  CLOCK  RATE'»/»'  2=1  MHZ ' r / r '  4=100  KHZ' r/f 
6=10  KHZ ' f / » '  8=1  KHZ ' » / » '  10=100  HZ'f/»'  14=60  HZ'/) 


READ ( 5  *  5020 )  IMCR 
WRITE  < 7  » 5070  ) 

FORMAT ('  '  »  '  INPUT  NUMBER  OF  CLOCK  PULSES  BETWEEN  SAMPLES  ) 

READ ( 5  r  5020  >  JSR 


WRITE  ( 7 1 5080  )  ^  vim 

FORMAT (  '  ' t 'CORRELATION  COEFF.  CALCULATED?  0=N0* 


1 =YES ' ) 


READ  <  5 » 5020 )  ICALC 


SET  UP  LPS  CLOCK  STATUS  REGISTER  PARAMETERS 


JCR= ( IMCR+257 ) 

CALCULATE  ACQUISITION  FREQUENCY 

IF  (IMCR.EQ.2)  FREQ= 1000000 . 0 /FLOAT ( JSR ) 
IF  (IMCR.EQ.4)  FREQ= 1 00000  » O/FLO AT ( JSR ) 
IF  (IMCR.EQ.6)  FREQ= 10000 , O/FLO AT  <  JSR ) 

IF  (IMCR.EQ.8)  FREQ=1000 , O/FLO AT ( JSR  > 

IF  (IMCR. EQ. 10)  FREQ= 100 . O/FLO AT ( JSR ) 

IF  (IMCR. EQ. 14)  FREQ=60 .0/FL0AT( JSR) 
IFREQ=IFIX ( FREQ ) 

ACQUIRE  DATA  WITH  SUBROUTINE  AQUIRE 
WRITE(7f 5090) 

FORMAT ( '  CR  TO  START  ACQUISITION') 

READ ( 5  ? 5020 )  I 


- 


204 


CALL  AOUIRE 

C 

C  CALCULATION  OF  MEANS  &  CORRELATION  COEFFICIENT 

C 

RNP=FLOAT ( JNP ) 

R=0 . 33333E+33 

SUMA=0 . 0 

SUMB=0.0 

SUMA2=0 . 0 

SUMB2=0 . 0 

SUMAB=0 . 0 

DO  50  1=1 » JNP 

SUMA=SUMA+FLOAT ( I A ( I ) ) 

SUMB=SUMB+FLOAT ( I B ( I > ) 

50  CONTINUE 

IF  (ICALC.EO.O)  GO  TO  80 
DO  60  1=1 f JNP 

SUMA2=SUMA2+FL0AT  { I A  ( I ) >**2 
SUMB2=SUMP2+FL0AT ( IB ( I ) >**2 
SUMAB=SUMAB+FLOAT( IAf I ) ) AFLOAT < IB ( I ) ) 

60  CONTINUE 

SA=DSQRT  (  <SUMA2-(SUMA**2/PNP>  )  /  <.  RNF'-l  ,0)  ) 

SB  =  DSQRT  <  ( SUMB2-  ( SUMBf  #2/RNF' )  >/(RNP-l  .0)  ) 

SAB= ( SUMAB- ( (SUMA*SUMB)/RNP) >/<RNP-l .0) 

SASB=SA*SB 

IF  (SASB.EQ.O.O)  GO  TO  70 
R=SAB/SASB 

70  IF  (SASB.NE.0,0)  GO  TO  80 

R=0.33333E+33 
80  SUMA=SUMA/RNP 

SUMB=SUMB/RNP 
SUMAB=SUMA/SUMB 
C 

C  OUTPUT  DATA  TO  TERMINAL 

C 

WRITE ( 7 » 6000 ) 

6000  FORMAT < '  MEAN  A  MEAN  B  MEAN  A/B 

WRITE (7 » 6020) SUMA-SUMB- SUMAB ?R 
6020  FORMAT ( '  ' »E10 . 4 ? 3X » E10 . 4 , 3X t E10 . 4 » 3X • El 1 . 4 ) 

90  WRITE ( 7 1 6040 ) 

6040  FORMAT ('  '»///'  STORE  DATA  =1 »  RERUN  =2>  EXIT  =3') 

READ ( 5  f  5020 )  J 
GO  TO  ( 10  >  20  r 30 )  J 
C 

C  FORMS  DATA  FILE  ON  DISC 

C 

10  WRITE ( 7  r 6060 ) 

6060  FORMAT  < '  ENTER  DE  *  FILE . DAT  FOR  STORAGE'//) 

CALL  ASSIGN (2t 'DKJFTN2.DAT' r-1 t ' NEW ' » ' NC ' f 1 ) 

DO  200  1=1 > JNP 
A(I)=FLOAT(IA(I) ) 

B<I)=FLOAT(IB(I) ) 

200  CONTINUE 

WRITE<2) (A(I)  »I  =  1 fl024) 

WRITE(2)<B(I)»I=1»1024> 

ENDFILE  2 
CALL  CLOSE ( 2 ) 

GO  TO  90 
30  STOP 

END 

* 
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